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Chapter I 
INTRODUCTION 

« 

Recent investigation by Croswell and Chatterjee (1972) has shown that 
a small homogeneous dielectric sphere -may be used to focus the radiation from 
an open-ended waveguide into a relatively directive beam, and hence would be 
useful as an antenna. 

Historically, large spherical lenses have been of interest for microwave 
antennas because their symmetry permits the formation of beams in any direction 
by displacement of source alone. Perhaps the best known lens of this type was 
introduced by Luneberg (1944). It has a radially varying dielectric constant and 
will focus rays emanating from a point source on the surface into a collimated 
beam of parallel rays at the opposite side of the lens. The Luneberg lens, how- 
ever, is difficult to fabricate in practice. A survey of other inhomogeneous 
spherical lenses is given by Johnson (1964). 

The imaging properties of the homogeneous dielectric sphere has been 
investigated by Bekefi and Farnell (1956), using the methods of geometrical 
optics and the diffraction theory of optical aberrations. Their results suggested 
that spheres of small refractive index, less than 30 wavelengths in diameter, 
may be used to produce good diffraction images of point sources. 

Assaly (1958) measured patterns for heptane (c * 1. 94) filled spherical 
shells 13. 3 and 33. 3 wavelengths in diameter In this experiment, the lens was 
supported by a conducting screen and illuminated by a waveguide displaced sev- 
eral wavelengths from the sphere. Data for a 12. 4-wavelength diameter sphere 
(e^ = 2.34) has also been presented by Cheston and Luoma (1963). Their results 
confirmed that the homogeneous dielectric sphere may be used to produce 
desirable radiation patterns despite the presence of a large amount of spherical 
aberration. 

Using a somewhat different arrangement, Meinke (1970) measured patterns 
and gain for a monopole formed on a small hemisphere (less than 1-1/4 wave- 
lengths in diameter) over a ground plane. In his work, using dielectric constants 
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between 2.2 and 11.5, Meinke noted that the gain tended to oscillate 

* s 

around that which could be obtained using a uniformly excited circular 
aperture, due to resonances inside the hemisphere. 

Finally, recent measurements by Croswell and Chatterjee, some of which 
have been published (1972), have shown that waveguide excited plexiglas 
(e r = 2.57) spheres 2-4 wavelengths in diameter exhibit properties 
desirable for feed antennas used to illuminate paraboloidal reflectors. 

Although considerable attention has been given to the study and 
measurement of scattering from dielectric spheres (see, for example, Inada 
and Plonus, 1970; Barriek, 1970 ), no rigorous treatment using electromagnetic 
theory has been given to the spherical homogeneous dielectric lens as an 
antenna element. 

In this presentation, the radiation from a source in the presence of a 
dielectric sphere is treated as a boundary- value problem, and the solution 
is obtained using the dyadic Green's function technique. Particular emphasis 
is placed on the Huygens' source which is used to model the waveguide excitation 
of spheres. 

Throughout this work the source, together with the dielectric sphere, is 
considered as an antenna, and hence the properties such as radiation resistance, 
directivity, and antenna efficiency are studied as well as the radiation pattern 



Chapter II 

DYADIC GREEN'S FUNCTIONS PERTAINING TO A SOURCE 
IN THE PRESENCE OF A DIELECTRIC SPHERE 

2. 1. Introduction 

The effect of a dielectric sphere on the radiation from an antenna may 
be investigated as a boundary-value problem under the assumption that the cur- 
rent distribution is a known function of position. The solution to this problem 
may be obtained by the technique of the dyadic Green’s functions. 

In this chapter, the boundary conditions are first presented. Then, 
starting with the free-space dyadic Green's function in the spherical system, the 
functions of the third kind, G^, are formulated using the method of scattering 
superposition. Subsequently the dyadic Green's functions of the fourth kind, G^, 
are derived from G^ using the principal of duality. In the next chapter, G^ and 
G^ will be used to obtain expressions for the fields due to electric and magnetic 
current sources, respectively. 

Although the functions presented here form the groundwork for the 
derivation of the fields in Chapter El, they should also find application to most 
problems dealing with dielectric spheres. 

2. 2. The Dyadic Green's Functions Pertaining to an Electric Current Source 
in the Presence of a Dielectric Sphere 

Figure 2-1 shows the geometry of a dielectric sphere with an externally 
located electric current source. The regions exterior and interior to the sphere 
are characterized by different constitutive parameters and, therefore, different 
propagation constants k^ and k^, where 

k l i € i (Region I, R > a) 

k 2 = U / M 2 e 2 (Region n, R < a). 

Here and e_ denote, respectively, the permeability and permittivity of the 
sphere, which may be complex. For most applications. Region I will be free 
space, and hence and = e^. 
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observation point 



FIG. 2-1: GEOMETRY FOR A CURRENT SOURCE IN THE 
PRESENCE OF A DIELECTRIC SPHERE. 


With the electric current source J(R') located in Region I, the electric 
fields in Regions I and n may be expressed by 


E(R) * iu/i 1 



(r|r>)- JCROdV* 


Region I 


( 2 . 1 ) 


E(R) *iwp 1 ^^G 2 2l) (R|R')' J(R')dV' Region n (2.2) 

where G„(R IR 1 ) is the dyadic Green's function of the third kind and the harmonic 
time factor e has been suppressed. The first and second numbers in the 
superscripts of the Green's functions in the above equations denote the regions 
of the observation and source points, respectively. Those functions are solutions 
to the vector wave equations 
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✓ xV*G 3 -k^G 3 =i6(R-R*) , R > a 

VxVxG 3 - - 0 , R < a . - (2.3) 

Besides satisfying the radiation condition at infinity. 


lim R[vKG 3 (R|R‘)-ik^ xG 3 (R|R'j] =0 , 

R —>00 

G 3 also satisfies the boundary conditions at the surface of the sphere 

nxG (r|r')| * n X G„(r|r')I 

|R=a_ |R=a + 

— fix VxGJr|r')| * — nx VxG„(r|r‘)I . (2.4) 

M 1 6 |R*a_ 14 |R=a + 

The simplest way to find G 3 (r|r') is to use the method of scattering 
superposition in which we let the fields consist of a sum of incident and scattered 
waves. With this method, the form of the scattered waves are first constructed 
using eigenfunctions which will satisfy the boundary conditions at the sphere sur- 
face as well as having the proper form at R * 0 and R = oo . The coefficients of 
the scattered waves are then determined by matching the boundary conditions on 
the surface of the sphere. 

Thus G 3 (r|r') is treated as consisting of two parts 

G^ n> (R|R’) * G 0 (r|r , )+G 3s (r|r>) , R>a (2.5) 

G 3 21) (R |r') * G^^Ir') , R < a . (2.6) 


Here Gq(r | R') denotes the free space Green's function pertaining to an infinite 

region with the same constitutive parameters as Region I, and G (R |r') repre- 

* =b/21) 

sents the part of the wave scattered from the sphere. G 3g Cr|r') is used to 
denote the Green's function internal to the sphere since only a scattered field 
will be present in this region. 

The free space Green's function for this problem is given by Tai (1971) 



6 


as: 


CD n 
ik, V-, r-, 


n“l m=0 


COM' (kJ+Ng 1 ^ (kjNi CO , R > R* 
“mn 1 ®mn 1 1 gmn 1 ’ 


where 

and 


M p k,) + N e (kW 1 * CO , R < R’ 

gmn 1 ^mn 1 e mn 1 ®mn 1 
o o o o 


r 

U, 


m * 0 

m fO 


M (k) * j 


o mn 


n <kR) [? 


(2.7) 


9P (cos 0) 

m ^m, M sm n cos 

— : - r P (cos0) mpO - — . mi 

smd n cos r 90 sm 


( 2 . 8 ) 


n q „(k) * 


n(n + 1) 


j (kR)P™(cos 0) ^ m^R + 


cos 


emn 


kR n n 


sm 


■ _L JL 

kR 9R 


[RJ„(kRj] 


r " 9P m (cos 0) 

n cos .* 


„m. 

P (cos 0) 
n 


90 


. m6d +m . _ 
sm r ^ sm0 


mi 


(2.9) 


The two sets of spherical vector wave functions, represented by M and 
N in (2.8) and (2. 9). are solutions to the homogeneous vector wave equation 

2 

V*VxF-kF*0 

as shown by Stratton (1941, p. 415), and also satisfy the symmetrical relations 

N e Oc)=^VxM p (k) 

o mn k , §mn 


M e (k) = TJ VX N (k) . 
®mn k emn 


(2. 10) 


In (2. 8) and (2. 9), j^(kR) is used to denote the spherical Bessel function 

of order n, and P m (cos 0) represents the associated Legendre function of degree 
n 
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n and order m. The expressions in (2.7) with superscript (1) are obtained from 

(2. 8) and (2. 9) by replacing the spherical Bessel functions with spherical Hankel 

functions of the first kind, n ; (kR). A prime is used on a function to denote that 

n 

it is defined with respect to the primed coordinate variables (R*, 6', <f>') pertain- 
ing to the R* vector. 

In constructing G^(R |r') from the wave functions we make the following 
observations: 

1. For the source located exterior to the sphere, the posterior parts 
of G^^CrIr 1 ) and G^^CrJr') must be the same as G (r(r') 
evaluated at R < R*. This is necessary for us to match the boun- 
dary conditions at the sphere surface. 

2. The propagation constants ki and kg must be used in the anterior 
parts of G^^'CrIr 1 ) and Gg g ^(R |r’)« respectively, because of the 
regions in which the observations points are located. 

3. The spherical Hankel functions must be used in the anterior parts 

w(l 1) — I — 

of G„ (R R’) because these functions represent outward 
3s 1 

traveling waves. 

4. Finally, since the fields are finite at the origin, the anterior parts 
of g!?^(r|r*) must be constructed using the spherical Bessel 

jS 

functions. 


In keeping with these requirements we let 


G 3s 1> ®I R,) 


°o n 

_i V y\o « i 2n + l (n-m)! 

4t r 0 ; n(n + l) (n + m)! 

n=l m=0 


[■ 


a^ 1} (k )M ,(1) (k )+bg 1J (k )N* (1) (k )~| 

§mn 6mn 1 emn 1 o mn o mn ^ § mn 1 -* 


( 2 . 11 ) 


and 
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oo n_ 

1 y” 1 to. x \ 2n + 1 fa-m)! 

4jt 0 ; n(n + l) (n+m)! 

n=l m-u 


[■ 


41“ (k2)< ( ^i) + d" n N <k 2 ) N ;« 

O o o o o o 


* 1 >] 


(2. 12) 


The scattering coefficients a, b, c and d in the above equations, in 
general, will be complex quantities and represent the magnitude and phase of the 
sphere's contribution to the total fields. The superscript (1) on these coeffi- 
cients denotes that they are defined for the source located in Region I, and 
should not be confused with the meaning of a superscript on a spherical wave 
function. 

Since the tangential components of the resulting E and H fields must 
be continuous across the sphere surface, we apply the boundary conditions (2. 4) 
which for our problem become 


RxG^CrIr') *RxG* 21) Cr|r') 
-7-rxvxg^CrIr') - — rxg^CrIr')., 


rR = a . 


(2. 13) 


1 


Solving the above equations for the coefficients, we find a^ =* a^, 

emn n 

K (i) , (i) (i) (i) m (i) ° 

b_ * b , c p * c and d Q * d , where 
emn n ^mn n emn n 

O o 0 


k 2 PiiW]' k l K'nM' . 


(1) 


Pc 


w 


a 


n 


k l p 2>n W . . . k 2 L P 2Wj 

5 *! 


h (pJ 
n 1 


k l r p 2 i n (p 2 ) ]' , , k 2 ElW]' . 


(1) “l 


w - 1 


2 P 1 


W 


k 2 KWl 


"2 P 1 


k i hWl 

Pc 


- 2 


h (p ) 
n 1 
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k i IfiWl . . . k i [ p i J n <p i > ]’ 

j ( n ) - " 

(i) “i p i " 1 p i p i 


W 


k i hVV 

P, P, 


PjJtpJ 


2 

“ 2 P 2 P 2 


h n (P 9 ) 

n 2 


,(1) M i 


k l [ P l h n (p l>]' . . , k ! P’lWT , . , 

— — i — 3„(p 1 )- — ■ L — : — ■*- hjp,) 


n l m, 


* 'r i # 

n 1 


n 


P lWj . 


W-JT 


1 L P 2 i n (p 2 ) J 


h n (p l ) 


(2. 14) 


Here, 


p i" k i a 
p 2 “ k 2 a 


K w l“ ^K W T' 


U), 


where z (p) is either j (p) or h (p). Also, h (p) denotes h (p), since only 
n n n n n J 

the spherical Hankel functions of the first kind will appear in this work. 

For a dielectric sphere in free space, 

^ 1 m ^2 m ^0 ’ 


6 1 * e o J 

V € 0 “ £ r ' 

where is the relative dielectric constant and may be complex. Equations 
(2. 14) for this case become 


.( 1 ) 


P 2 i „ (p 2 y ] 


p lW] 

' in (p 2> 

hW! 

’W* 

p 2 1 n <P 2 ) ! 

'V P 1 > 


d d> _ '/s [ p i h l x V]' W- ^7 EiWTW . 

e rf P l h n 1)( ' > l ) ]' J n <P 2 > - r p 2 J n fe> 2 ) ]’\ ( ' > l ) 

(2. 15) 

Substitution of equations (2. 1), (2. 11) and (2. 12) into (2. 5) and (2. 6) 
gives us the desired Green’s functions for the configuration of Fig. 2-1. They 
are 


Source in Region I, Observation Point in Region I 


^ (11) (r|r) * — - 


oo n 


/O r X 2n + l (n-m)! 

‘3 '“TV 4 jr ^ — i'^"*Vn(n + l) (n+m)! 
n=l m-u 

Mg X) (k ) Fm’ (k.J + a^M'^ (k )1 
®mn 1 L e m n 1 n e m n 1J 


+ N^ 

o mn 




gmn 1 


.. 00 n 

ik, \ 1 v 1 


(kJ+b^Ni^Ocj) 


gmn 


(11) mlR’> * — i > ) (? , X 2n + l (n-m)! 

3 ®I R) in ^J 2 - 6 o’ n(n + l) fa+m)! 

n=l m=0 

' {[“e mn ft 1 ) + % , «elfl< 1 >] ft" 0^ 

+ ["n (k )+b (1) N (1) (k )1 n' ( 1) (k )\ , 

L ®mn 1 n e 1 J emn 1 l 

o 0 mn 0 J 


R > R’ 


(2. 16) 


a < R < R' 


(2. 17) 
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Source in Region I, Observation Point in Region n 

., oo n 

G (21) (rIr') * — ^ ( 2 - 6 ) — - + -- . 

3 ' 4^ L - — i — • 0 n(n+l) (n+m)! 


n=l m=0 


Fc (1) M e (k )M» (1) (kJ+d (1) N (k_)N' (1) (kjl , 

I n mn 2 emn 1 n e mn 2 e mn 1 

u r\ n r\ n 


R < a . 

(2. 18) 


For the source located inside the dielectric sphere, again consider 

x(22) — I — 

the geometry of Fig. 2-1, but now with R' < a. We let (R|R') represent 
the free space Green's function in an infinite region with the same constitutive 
parameters as that of Region n. From (2. 7) we have 

oo n 

G (22, <rIr-> - -2 y y . > Jn±i- is^ssll . 

0 • 47r ' 0 n(n+l) (n + m)! 

n*! m=U 


r M^ 1} (kj M' e (kJ+N^ (kJN' (k ) , R>R' 
®mn 2 ®mn 2 emn 2 e mn 2 


^ (kJM' (l) <kJ+N a (k_)N' (1> (kj , R<R' . 


emn 2 emn 2 e mn 2 ®mn 2 
o o o o 


(2. 19) 


=( 22 ) — — 

Using the method of scattering superposition, G^ (R R') is treated as con- 
sisting of two parts, 

G^ti* J5') » Gq 22) (R |r') +G 3 (22) (R |r«) , R < a 


G< 12) (R |r ') - G‘f(R |r') 


R > a 


( 2 . 20 ) 


=( 22 ) 


Here G„ is used to represent the portion of G„ which is scattered from the 

dS =(12; 
inside surface of the dielectric sphere, and G is the transmitted field. We 

s (22) *(i2) 

note that G^ s must be finite at the origin, must satisfy the radiation con- 

dition at infinity, and the boundary conditions must be met at the sphere sur- 


face . Thus , 
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P (12) (rIrM = (2-6 ) --- iSLZJ93>i 

3s I 4?r , ^ 1 0 n(n + l) (n+ m)! 
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[■ 


e 2> Mg 1} (k )M^ m (k )+b^ 2) Ng 1} (k )N' (kj] , 
®mn gmn 1 0 m n 2 gmn gmn 1 gmn & J 


R > a 


( 2 . 21 ) 


and 


Gf/Vk) 


ik 


oo n 


2 

4ir 



n=l m=0 


(n r \ 2 n + 1 (n-ra)! 

0 n(n + l) (n + m)! 


r4 2) M e (k )M' (k ) + d^ N e (k ) N' (k )"] , 
L ®mn ®mn 2 e mn 2 e mn « mn 2 e mn 2 J 


R < a 


( 2 . 22 ) 


where the superscript (2) on the scattering coefficients denote they pertain 
to the source located in Region II. 

The coefficients are now found by applying the boundary conditions at 
the sphere surface, namely 


A 3B 

R*G 


^ 12) (r|r') = RxG 


? 2) (r|r’) 




1 /\ = 
— R x VXG 


( 12 ) 

3 




R * a . 


We again find 


,® 

o mn 



b® 

emn 

o 



(2) 

c e 

o mn 


( 2 ) 

c 


and d 


( 2 ) 


§mn 



where 


(2) hWT W - r p 2 h n (p 2? i . (p 2 ) 

3 . 

[f2 J n <p 2>]' h n (p l > - ^ D°l ll n <p l ) ] W 



[ p 2 h n (p 2 ) ]'in (,i 2 > 


V2 


t?2 i n <P 2 > j' W ' ^ [ p l h n (p lD' W 
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For the case where n * /u , these equations become 

1 £ 


(2. 23) 



Substitution of equations (2. 19), (2.21) and (2.22) into (2.20) yields 
the desired Green*s functions for the internally located electric current source 
as follows: 
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Source in Region n, Observation Point in Region I 

.. oo n 
lk 




(2 - 6 ) 211 + 1 fa-m)? 

4tt tri* 0 n fe + 1 ) (n+m)! 

n=l m-u 


• R>a < 2 ' 25 > 


Source in Region IT, Observation Point in Region n 


oo n 






n=l m*0 


Mp 1} (k„)+c (2) M e (k )*] M ' (kj 
e mn 2 n e mn 2J e mn 2 
o o o 


n*2>} • ^ R>E ' 
° o O J 


( 2 . 26 ) 


k „ 00 , n 

5< 22) (5|5') - C-V X *\ ) “ • 


n*l m=0 






M ~ ^ > L M emn <k 2 )+C n M e 


mn 


(k ? )+d l 2>N 4 ^9>lr ’ 

o mn 2 L* gmn 2 n gmn 2 -* J 


R < R' . 

(2. 27) 


2 . 3 Duality Principle Applied to and 

When dealing with magnetic current sources, it will be convenient to use 
the dyadic Green's functions of the fourth kind, G^, to obtain the magnetic field. 
The symmetry which exists in Maxwell's equations provides a simple means 
of obtaining these functions from those of the third kind which were derived in 
Section 2. 2. 
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Consider the following sets of Maxwell's equations: 

Vx E * itquH 

V X H * -iucE + J 

V-J«iwp (2.28) 

VX E 1 * iupH* - J 

m 

Vx H' * -iioeE' 

V * J * iiop (2. 29) 

m m 

where and p^ are the magnetic current and charge densities. Here, the 
primed and unprimed fields are due to electric and magnetic sources respectively. 

As shown in detail by Papas (1965), a duality exists in the above equation 
sets. This principle, however, may be extended to include the dyadic Green's 
functions of the third and fourth kinds by considering the inhomogeneous vector 
wave equations derived from (2. 28) and (2.29) and their solutions. Thus, 


and 


_ 2 - - 
VxVxE-k E * iupJ 

Vx VXH-k 2 H * koeJ 



(2. 30) 


(2.31) 


or 


E =iwp^G 3 (R|R’)'KdS’ 


H' = iue UG.CR R')*K dS* 
4 m 


(2. 32) 
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where the surface electric and magnetic current densities are given by 

K * -n x E' 
m 

K * n X H . (2.33) 

To change from a system of fields excited by electric sources to one 
excited by magnetic sources, or vice versa, we see from equations (2. 28) to 
(2. 33), it is only necessary to replace the quantities by their duals as sum- 
marized in Table 2-1. 


Table 2-1; Dual Quantities 


E 

H 

S 3 

J 

K 

P 

p 

€ 

k 










H* 

-E' 

°4 

J 

m 

K 

m 

P m 

m 

e 

P 

k 


It also should be noted that besides satisfying the radiation condition, 
also satisfies the boundary conditions at the surface of the sphere. 


n xG 4 (R|R') 


R=a 


n XG 4 (R|R') 


R-a. 


VxG 4 (r|r’)| 


|R=a 


» — n x Vx G (R |r') 
e 2 4 


(2. 34) 


R=a . 


which are the dual of (2. 4). 

2.4. The Dyadic Green's Functions Pertaining to a Magnetic Source in the 
Presence of a Dielectric Sphere. 

Since the propagation constant k is its own dual, it can be seen from 

(2. 8) and (2. 9) that the spherical wave functions will remain invariant in the 

transformation from G^ to G^. The only difference, therefore, will be manifest 

in the scattering coefficients a , b , c , d , and once these coefficients have 

n n n n 

been determined, it is possible to list the various Green's functions of the fourth 
kind. 

35 

To obtain the coefficients for G 4 from those previously derived for 
G 3 we let 
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a* (p, e) * a (e,p) 
n □ 

b’(p,e) *b (e,p) 
n n 

c^(iu,e) * c n (e,A<) 

d' <p, e) * d (e,p) 
n n 

3C 

where the primed coefficients are used in G, and a , b , c > d are defined by 

4 n n n n 

(2. 14) or (2. 23) depending on the source location. Equations (2. 15) and (2. 24) 
were not used for the unprimed coefficients in this case because p^ cannot be 
set equal to until after the duality principle has been applied. 

From (2. 14) for the external source we find 
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b< 0 '- 
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P 2 j n <P 2 ) ] 


(p l , - £ rE 

“lW]' 


£ rl P l h n <f 

( p 2 ) n <P 2 ) . 

V] 

l^s 

'W- 

W-"r| 

t p 2 ) n (p 2 ) . 

flW]' 

|' h n (p l ) 

W 

^rpl h n (p l[ 

'W* 

P 2W" 

fw 



e 

r 


’j (p )-e 
J n 1 r 


'W 

e 

r 

fiVV 


, p 2^n (p 2 ) _ 

*1 




n 


k l k 2 C r [ P l h n (p l > T W ' k l k 2 £ r[ p l i n (p l ) ]' W 
k 2[ p lWl’ W - k l £ r f 2 i n (p 2 ) ] ' W 

(2. 35) 


where we have used the fact that ^ 1 /p 1 = kg/P 2 * aud let P r = For the 

case when p^ = Pg, equations (2. 35) may be further simplified; it is found that 


(U* - h ( 1) 
a = b 
n n 


. UV - J 1 * 

b — a 
n n 
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(2. 36) 


Proceeding in a similar manner for the internally located magnetic 
source, it is found from {2. 23) 


^2^2*7 h n (p 2 ) " 


! h n (p 2 ) ]'^n (p 2 ) 

P 2 i n <P 2 ) ]' h n <P l ) " 

^1 

?lWj j n <P 2 ) 


(2)' f^n^V] * V p 2* ~ P 2V P 2Q tj n (p 2 ) 

U 

[ p l h n (p l0' i n (p 2 ) 



j(2)' 

a 

n 


e r[ p l h n (p l > ]' h n (p 2 ) - 

P 2 h n (p 2 ) ] 

'W 

t P 2>n ( ' , 2U , W- £ r : 

P l h n (p l>] 

'W 


€ r[ P 2 i n <P 2 ) ]'V P l ) ‘ (^) [ p l h n (p l , ]' i a (p 2 ) 


(2. 37) 


By letting it is apparent from the above equations that 


( 2 )' 

a 

n 



b< 2) 



1 ( 2 ) 



n 


,( 2 )» ( 2 ) 

d * c 
n n 


(2. 38) 


where the unprimed coefficients are given by (2. 24). 
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It is now a simple matter to list the dyadic Green's functions of the 
fourth land by making the appropriate substitution of the scattering coefficients 
in equations (2. 16) to (2. 18) and (2. 25) to (2. 27). These functions are listed 
below for the various regions of source and observation points. 

Source in Region I, Observation Point in Region I 


oo n 


(m^OcjTm' (k )+b (1) M' (1) (k.il 
) e mn 1 L e m n 1 n e m n 1 J 
o o o 

+ (k ) [Fji (kj +a (1) N4 (1) (k.)"| 

§mn 1 L§ m n l n ®mn 1J 


R>R» 


ik 


oo n 


4 u) (r|r’>“ ' 2 - 6 «> 


n=l m=0 


2n + l (n-m)? 
O' n(n + l) (n+m)! 


[m (kJ+bJ^M^ (k )]m« ( 1) (kj 
L ®mn 1 n £mn lJ e mn 1 


+[n g (k ) +a^ 1) Ng 1) (k )]N^ (1) (k : 

L o ma 1 n §mn 1 J £ mn 1 
Source in Region I. Observation Point in Region II 


> , a < R < R' (2.39) 


oo n 


g (2i) (rIr') * y 1 (2-6 ) 2n+i . 

4 ' 47 r 0 n(n + l) (n+m)! 

n*l m=0 

• [d^Mg (k )M' (1) (k ) +c^N (k > N' (1) (kj\, R < 
L n e m n 2 e mn l n e mn i e mn 1J 


(2.40) 
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Source in Region It, Observation Point in Region I 


ik. 


oo n 


r (12 ) ( r|r«) * — /> y (2 6) 2n + 1 * n ~ m)! . 

G 4 ®I R ' in S—f L-, < 2 - 6 0 > n(n + 1) fe+m)! 
n=l m=0 

• |V 2) Mg X> (k )Mg (kJ + a^Ng (k )N* (k )~| , R>a 
L n o mn 1 2 n omn 1 gmn 2J 


(2.41) 


Source in Region n. Observation Point in Region H 
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(k 9 ) +cf } N (k„)l N' (kjl , a > R > R 1 
U omn 2 n gmn 2J e mn z 


oo n 

g< 22) <r|5') - ^ Z £ < 2 - v “ 


n=l m=0 


4 


ML (1) OO +d <2) i\i; (k 2 ) 


o ma ^ 


n e 


(kjfi 

mn ^ 

KJ 

+n (k )fe (1) (k )+c (2) n; (k )1 
gmn 2 L gmn 2 n gmn 2 J 


>, R<R' 


(2. 42) 


where we have used the fact that for k^ * k^/ \j~^‘ 

This completes the derivation of the dyadic Green's functions, pertaining 
to a homogeneous dielectric sphere, which will be used in subsequent chapters. 



Chapter HI 


ELECTRIC FIELDS DUE TO SIMPLE POINT SOURCES 
IN THE PRESENCE OF A DIELECTRIC SPHERE 

3. 1. Short Horizontal Electric Dipole 

Consider an infinitesimal horizontal x-directed electric dipole with 
current moment C g located at R' = b, 0' =0, {>' =0 as shown in Fig. 3-1. 


z 



FIG. 3-1: HORIZONTAL ELECTRIC DIPOLE IN THE 
PRESENCE OF A DIELECTRIC SPHERE. 


We let the electric current density be represented by 

J (R') = C . (3.1 

6 b sin0' 

To find the electric fields in the two Regions, (3. 1) is substituted into 
(2. 1) and (2. 2). Because of the delta functions within the integrand, the follow- 
ing is obtained 
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E(R) *iuu C G< U) (r|r') 4 

1 6 o 


E(R) *iu^C g1 21) (r|r')*'x 
led 


R ,as b 

e*=o 

^’*0 ‘ 


R’*b 

e‘*o 

^'*0 


Region I 


Region n. 


(3.2) 


The expression G^CrIr 1 ) must now be examined. First, the 
associated Legendre functions in (2. 8) and (2. 9) are evaluated at 0 = 0. 
It is known (Harrington, 1961) that at 0*0, i. e. , cos 0 * 1, 


P m (l) 


1, m * 0 
0, m ^ 0 


Also, using de l'Hopital’s Rule for the limit. 


p m (cos 0) 
n 

sin0 


0*0 


f n(n + 1) 

1 2 
l 0 , 


m * 1 
m ^ 1 


These last two equations show that for the point source located at O’ * 0, only 
those terms in the series expansion corresponding to m * 1 will contribute to 
the fields. 

Since in the spherical coordinate system 

x * sin0cos^R + cos0cos^0-sin^ , 


it is apparent that 


M (k) -x 
emn 


j £=0 


N (k) 
omni 


j ^*0 


and hence 


M, 


e <k) -x 
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M . (k) 
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n(n + 1) 


j n (kb) 


(b, 0,0) 


l(b, 0,0) 
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N, 


e 


(k) 


•SI 


mn 


"ela 00 


Kb, 0,0) 


n(n + 1) 
2 


[kbj n (kb)]* 


kb 


(b.0,0) 


(3.3) 


where the M and N functions are given by (2. 8) and (2. 9). 

Substitution of the dyadic Green’s functions given by (2. 16) to (2. 18) 
into (3. 2), and using the results (3. 3), yields the following expressions for the 
electric fields: 


Source in Region I, Observation Point in Region I 


- - _ jge f 2n + l 

E e®* ~ 4tr •'— ! n(n + l) 
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a < R < b 


(3.5) 


Source in Region I, Observation Point in Region n 


E (R) = 
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(3.6) 
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Here, = k^b and the subscript e is used to denote that the field is due to an 
electric current source 

For the current source located inside the dielectric sphere, consider 
the geometry of Fig. 3-1, but with b<a. The electric fields will then be given 
by 


E(R) =iup 2 \\\Gg 12 (R[R')-J(R’)dV' 


or 


E(R) = ii^ 2 \\\Gg 22 |R|R') • J(R')dV' 


(3.7) 


depending on the location of the observation point. 

Letting p = p 2 and using the dyadic Green's functions given by (2. 25) 
to (2. 27) together with (3.3), the following expressions are obtained. 

Source in Region II, Observation Point in Region I 
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Source in Region n, Observation Point in Region II 
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R < b (3. 10) 


where p „ = k„b. 

4 2 

3. 2 Short Horizontal Magnetic Dipole 

In this section, the radiation from an infinitesimal magnetic dipole located 

on the +z axis is considered. The source is pointed in the -y direction and has 

a magnetic current moment C . 

m 

The magnetic field is given in terms of the dyadic Green's functions of 
the fourth kind for the appropriate regions of source and observation points as 
follows : 


H(R) = iue 1 




J (R')dV' 


H(R) = iue 



5< 21) <r|r'W _(R')dV' 
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H(R) - iue 



p^ 12) (R R') • J (R')dV' 
4 m 



(3. 11) 


where J m is the magnetic current density, and for this problem is expressed by 
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J (R*) = -C 
m m 


6(R*-b)6(0')6(^') a 

2 ■ y 

b sin0' 


(3. 12) 


In the spherical coordinate system, 

y = sin0' sin^'R' + cos 0' sin^ , 0 1 +cos 


(3. 13) 


On performing the integration of (3. 11), one finds the magnetic field will be 
given by 


H(R) - -itoeG^CR |r') ‘ y 


(3. 14) 


From (2. 8) and (2. 9), it is seen 
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(3. 15) 


where we have used the relationships for the associated Legendre functions 
derived in Section 3. 1. 

The magnetic fields due to a magnetic point source could now be listed, 
but it would be more convenient to find the electric fields. In a source-free 

region, _ 

E = — V X H. 

U€ 


V *M 0 (k) = kN e (k) 


?mn 


mn 


Also, since 
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and 


VXN (k) - kM. (k) , 


emn 


;mn 


we simply need to multiply each term by the factor ik/ue appropriate to the 
region of observation, and interchange the roles of the resulting M and N 
functions to obtain the electric fields from the magnetic field expressions. 

Using the above result and the Green’s functions given by (2. 39) to (2. 42), 
we find the electric fields due to an infinitesimal directed magnetic dipole 
located at (R 1 =0, 0' = 0, <f > ' = 0) may be expressed as follows. 

Source in Region I, Observation Point in Region I 
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Source in Region I, Observation Point in Region II 
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Source in Region It, Observation Point in Region I 
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Source in Region n, Observation Point in Region II 
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Here the relation ri~ ^p/e is used and again = k^b and p^= kgb. 



Chapter IV 


ELECTRIC FIELDS DUE TO A HUYGENS’ SOURCE 
IN THE PRESENCE OF A DIELECTRIC SPHERE 


4. 1 Introduction 

The radiation from a homogeneous dielectric sphere, a few wavelengths 
in diameter, placed over the aperture of an open-ended waveguide is of consider- 
able practical interest since it has been shown experimentally (Croswell and 
Chatterjee, 1972) that this device exhibits properties desirable for a feed antenna. 

Because it would be exceedingly difficult to solve this problem exactly, 
we will neglect the scattering from the waveguide and model the exciting fields 
using a superposition of electric and magnetic dipole currents forming an 
approximate Huygens’ source. 

In this chapter, the electric and magnetic current relationship necessary 
to construct a Huygens’ source in free space is first derived, and then the radi- 
ation from it is examined. Subsequently, the electric fields are derived for an 
approximate Huygens’ source in the presence of a dielectric sphere, and it is 
shown that the same results are obtained if the source is allowed to approach 
the sphere surface from either the inside or outside of the sphere. Finally, the 
far-zone fields are derived and the results examined for a very small sphere. 


4.2 Huygens' Source in Free Space 

Consider a small portion of a uniform x-polarized plane wave in free 
space, incident normally on the z = 0 plane as shown in Fig. 4-1. The electric 



FIG. 4-1: PLANE WAVE GEOMETRY 
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and magnetic fields across this plane are related by 


E 

x 



(4. 1) 


When the fields to the right of z = 0 are of interest, the incident fields may be 
replaced by fictitious electric and magnetic current sheets. 


K = nXH 
e 

K = -nX E . 
m 

For this problem, n = z and thus 

K = -H x 

e • y 

= (4 ' 2) 

Allowing the portion of the plane wave under consideration to become 
small, the incident fields may be replaced by the equivalent electric and mag- 
netic dipole current moments C x and C y, where 

e m 


C = 
e 




(4, 3) 


Substituting (4. 2) into (4. 3) we find C * n 0 C . This double source corres- 

m 0 e 

ponding to the ratio of electric and magnetic fields of (4, 1) is called the "Huygens’ 
source" because Huygens introduced the concept of considering a wavefront as 
a system of secondary sources. 

Let us now derive the far-zone fields radiated from a pair of crossed 
point electric and magentic dipoles as shown in Fig. 4-2, and then particularize 
to the special case of the Huygens’ source. 

For the x directed electric dipole,the far-zone electric field is given 
by 
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FIG. 4-2: COORDINATE SYSTEM FOR HUYGENS’ SOURCE 


- e ikR 

E = - C ^z—sin 0'0' (4.4) 

e 4w ■ e R 

where C is the electric current moment of the source and 0’ is measured with 
e 

respect to the x-axis. From duality. Table 2-1, it is found that the magnetic field 
radiated by a -z directed magnetic dipole of current moment C^ is given by 



iue n C ikR 

0 m 6 sin0 M 0" , 


47 r 


R 


(4. 5) 


where, in this case, Q" is measured with respect to the y-axis. 

A — — 

In the far-zone R X E = rjgH, and hence from (4, 5) 


-imue C ikR A 

E = — sin0"p" . 

m 47T R 


(4. 6) 
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Using the coordinate transformation 

sinS'd' =-cos0cos^0+sin^$ - 

sin0 , '0” = cos <f>6 -cos 0sin^£ , 

where 0 is measured from the z-axis as is the usual case for the spherical 
coordinate system, (4. 4) and (4. 6) become 

iwAf C ikR A 

g __ — _ — - e (—cos 6 cos id + sin^ 2) (4.7) 

e 47T R r r r 


__ iwr? i e C ikR A 

E — — (cos id - cos 0 sin J>7>) . (4.8) 

m 47T R rv 

Letting 

C =n C , (4.9) 

m e 

where q = rjg //V (4.7) and (4.8) may be combined. Thus, 


E = E +E = - 
Hem 


iu/u n C ikR 
0 e e 


47 r 


R 


-cos ^cos)6e+ ^1 - sin^ 


(4. 10) 


For now the factor is simply a constant used in the ratio of the electric and 
magnetic current sources. 

We let e *= 1 to construct the Huygens' source in free space, and hence 
(4. 10) becomes 


__ iw/^C ikR A 

E H = ~ — 4ir 6 (cos^+sin^d-cosfl) . (4.11) 

It is apparent from this equation that the Huygens ' source will radiate a linearly 

A 

polarized cardioid-shaped pattern with a field maximum in the -z direction and 
a null in the +2 direction . 

It is because of this directional characteristic together with its 
simplicity that the Huygens' source has been chosen to model the waveguide 
aperture excitation of a dielectric sphere. 
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4. 3 Huygens* Source in the Presence of a Dielectric Sphere 

We can only approximate the Huygens' source in the presence of the 
dielectric sphere because the complex wave impedance r) is not known. 

Although, in general, r) would be a function of the sphere diameter and 
constitutional parameters as well as the source location, we will let 



b > a 
b < a 


(4. 12) 
(4. 13) 


where ^ 2 = \j 4*2' the re & ions are shown in Fig. 


4-3. 



FIG. 4-3: HUYGENS' SOURCE IN THE PRESENCE OF A 
DIELECTRIC SPHERE. 


Equation (4. 12) was chosen because it would be valid if b»a, since for 
this case the sphere would have little effect on the source. Alternately, (4. 13) 
would be appropriate for the source located within, or on the surface of, a large 
lossy dielectric sphere. Also, as will be shown in subsequent sections, (4. 13) 
yields a good approximation to the Huygens’ source on the surface of small 
spheres. 
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To obtain the electric fields due to the source located external to the 

sphere, let C = r\.C and combine (3.4) to (3. 6) with (3. 16) to (3. 18). A 
m 1 e 

similar procedure is followed using equations (3.8) to (3. 10), (3. 19) to (3. 21), 
and letting = r^C^. resu ^ing equations for the electric fields in the 
various regions are given as follows. 

Source in Region I, Observation Point in Region I 


-k.ujuC 


V 0- 


47 r 


1”" e y -1 2n + l 

+—i n(n + l) 
n=l 


E h (R> = 


+ 


h (p )-i 
n o 


[ P 3 h n (>, 3 ) T 


)M (1) (k.) + (B (1) +C 

oln I 

n i 

UpC <r°° i 

e \ 

2n + l 

47T 

n(n + 1) 

n=l 

r p 3 h n <P 3 ) ]’l 

IPs 

P 3 . 

[ L oln 


R > b (4. 14) 


w (k 

n oln 


+ih(Po) 
n o 


N (k )+b^N^ (k 
e In 1 n e In 


a 

,>]V 


Source in Region I, Observation Point in Region IT 

oo 


E h ®> 


-k,unC 
1 e 


V\ 2n + l 

n(n + l) 


47T *■ — i. n(n + l) 

n=i 


{[■ 


+ 


W- i p , 3 h n Cp 3 ) l '] c n >M o lA* + 

hW]’ , , , ; 

J= — +ih (pj 


n r 3' 


d (1) N (k ) 
n eln 2 


, 0 < R < a 


a < R < b 
(4. 15) 


(4. 16) 
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Source in Region It, Observation Point in Region I 


e h®> 


-k.upC 
1 e 


oo 

E 


2n + 1 


47r i n(n + 1) 
n=l 


((A (2) + D (2) )M (1 1 ) (k ) + (B (2) +C (2 ^)N^ (k )V R >a 

\ n n oln 1 n n eln 1 J 


Source in Region n. Observation Point in Region II 


oo 


2n+ 1 


_ _ - k 2 ^ c e 

E H <R) 4 v n(n + l) 

n=l 
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fn , . . hwT 

^ w - 1 p . ■ 


hW! , .. 


* l W 


(k ) +c®^M . (k ) 
oln 2 n oln 2 


N eln (k 2 )+d I\l» l V 
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a > R > b 


2n + l 


_ _ - k 2 W/jC e 

E «®> - -IT- t-i n(n + 1) 
n=l 


oo 

E 


fr 


n 4 n n 4 


<2), (0 [>M!f _ ic (2) fcwT 

m A ' n n A 


M (k r 
oln 2 : 


+ 


o. 4)+ m® ,(p 4 ) 

p„ . n p„ n 4 n J n v 4 

4 4 


N eln®2> 


R < b 


A n ' “ 'A 1 +a n 1 Vs* 

(1) NnM' (1) [ P 3 h n (p 3 ) ]' 
B — rb 


(4. 17) 


(4. 18) 


(4. 19) 


where 
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and 


,Ml! u (d fowl 

n P 3 1 n p 3 



(4. 20) 


(4.21) 


The coefficients in (4. 14) and (4. 17) were combined into A , B , C , 

n n n 

since the primary interest here will be the electric field far from the antenna. 

A and B represent the field contribution from the electric current source and 
n n 

similarly C and D were derived from the magnetic source. Again, the 
n n 

superscript on the coefficient denotes the source region. 

These coefficients are now simplified for the case where the source is 
located on the surface of the sphere. The procedure will be to first let the exter- 
nal source approach the sphere surface. Subsequently, we let the internal source 
approach the surface, and the results are compared as a check on the validity 
of the expressions. 

Since C = n,C for the external source, (4. 14) must be modified for 
m 1 e 

the source on the surface where C = rj n C . Thus, (4, 14) becomes 

m 2 e 
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E h® 


-k.u/uC 
1 e 

47T 


z 

nr 1 


2n+ 1 
n(n+ 1) 






b*a 


(4. 22) 


where the relation = rj 2 ^ e r ^ as been used. 

To evaluate the coefficients, we first insert (2. 14) into (4. 20) and let 


p = p , i.e., k b = k a. The Wronskian relation from Appendix A, 

O J. J. J. 


) n (p) [pb n (p)] ' - h n (p) [pJ n (p)]' = I/P , 


(4. 23) 


is then applied to simplify the expressions. After some algebraic manipulation, 
it is found that for the source on the sphere surface 


A (1 > - - 
n 


B (1 > = - 
n 


c* 11 

n 

VS 

<i) 


i 

W 

P 1 [ P l h n (p l>]' 

i n <p 2 ) -[ p 2 i n (p 2 ) ]' h n (p l ) 

i 

p 2 i n (p 2 > ' 

2 c fp h (p 
P 2 r(_ 1 n 1 

) ]' i n <p 2 ) -[ p 2 j n (p 2 , ]' h a <p l > 

e 

r 

W 

• jm% m -■ «-' 1 - ^ ^ 7~ ~~ 


D 


n 


Vs 


p 2 SflWll'W ~ PPSIFW 
1 

p l p 2 [sW} W ' [ p 2 i n (p 2 ) ]' W ' 


(4. 24) 


In a similar manner but with p^ = p 2 we find from (4. 21), (2. 24) and 


(4. 23) that 
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A< 2 > = A (1 > 
n n 




B (2> = B (1) 
n n 


c (2) = c {1) /JT 

n n ' V r 


> 


p l =p 3 


P 2 =P 4 


d ( 2) = D (1) /JT 
n n ' ¥ r 


J 


(4. 25) 


In comparing (4. 22) with (4. 17), it is seen that (4. 25) confirms that the 
same results are obtained if the excitation is allowed to approach the surface 
from either the inside or outside of the sphere. 

4.4. Far-Zone Field 

In the far zone, the spherical Hankel functions of the first kind may be 
replaced by the leading term of their asymptotic expressions. Thus, 


where 


, , ikR 

h n (kR) ~ (-1)“ \ 


® h KH ~ <-‘> n 


ikR 

i 

"kR 


(4. 26) 


For large values of kR, the spherical wave functions become 


ikR 

M (kR) ~ (-i) — - m 

gmn ^ Qmn 


Ng 1} (kR) ^ (-i) n - 
gmn 


ikR 


, ■ n 

kR e mn 


(4. 27) 


m 


^mn 


f ~ _J22_ p 1 * 
sin 8 n 


x ap m (cos e) 

m, sin ,£ n cos 

(cos 9) mpd - — . mi 

cos r dd sin 


] 


(4. 28) 


and 



39 


3P m (cos0) P m (cos0) 

n cos /^ - n sm 

rr . mpo 4 . ; — mp p 

30 sm r + sin0 cos yy 


(4. 29) 


Using these results in (4. 14) and (4. 17) we obtain the following far-zone 
field expressions for the Huygens’ source in the presence of a dielectric sphere. 

Source in Region I, C =n,C 
— “ 1 m l l— e 


-up,C ikiR r — 1 _ , , 

v „ - 1 e e L > - 2n + l / ,vD 

E H (R) 4tt R f n(n + l) * * 

n=l 

■T-i(A (1) +D (1) )m +(B (1) +C (1) )n 1 

L n n oln n n elnj 


(4. 30) 


Source in Region n, C 331 n„C 
— m 2r~t 


-u/i.C ikiR T — ■ _ . , 

w>“ 

n=l 


{- i( A <2 > + D (2) )m , +(B® +0®>n , \ 
^ n n oln n n elnj 


(4.31) 


where the coefficients are defined by (4. 20) and (4. 21). 


In order to recover the far-zone field expressions due to the electric 

or magnetic source only, simply let C = D = 0 or A = B = 0 in the above 

n n n n 

equations . 

Now the source is placed on the surface of the dielectric sphere and the 

diameter is allowed to become very small. The results using C = rj,C and 

m 1 e 

C = C are then compared to those of an ideal Huygens’ source in free 
m 2 e 

space. 

For small arguments only the first term in the ascending series for the 
spherical Bessel functions is needed. It is found (Abramowitz and Stegun, 1966) 


1*3*5... (2n+l) 


h (p) '\x 1 * 3 • 5. . . (2n- l)(-i)p 
n 


(4. 32) 
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hence. 


K (P 0’ ~ 1.3-5^(2n + l) p “ 

{ph n (p)"3' 1*3*5... (2n - 1) * i • np + ^ 


and therefore. 


[ P lV P l ) ]’ i n (p 2 ) '' / 

| p 2 J n <P 2 ) ]' h n <p l ) ~ 


in -(n+1) n 
2n + l P 1 p 2 

-i(n+ 1) -(n+1) n 

2n + 1 P 1 p 2 * 


(4. 33) 


Using these relations in (4. 24), 


A (1 >~ 


n 


1* 3*5... (2n + l) 


R (l) ^ (2n+l)(n + l) (n— 1) 

n (e n+n + 1) ‘ 3 • 5. , . (2n + l) P 1 
r 

(rt i€ < 2n + 1 ) 

r> *v> — n n 

n (e n+n + 1) * 3 * 5. . . (2n + l) P 1 


. (n-1) 

D (D 

n 1*3*5. ..(2n + l) ’ 


(4. 34) 


and therefore in the limit as p — >0 


A< U > C (1) 
n n 


0 



n 



2 

€ +2 
r 



n 



2i 

3 


(4. 35) 


Inserting (4. 35) into (4. 30) and evaluating the m and n functions given by (4. 28) 
and (4. 29) we obtain the following expression: 
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E H (R>= 4* 


-*’*‘l C e e ik l R 


R 


ft'^ 0086 ) cos j>Q + (- "T77 - cos 0 


k€ r + 2 


0^sin$i/^ , 


(4. 36) 


1 3 1 

where the facts that P (cos 0)/ sin 0 = 1 and — P (cos0) = cos 0 have been used. 

1 da 1 

Had C = n 0 C been used in relating the electric and magnetic current sources, 
m 2 e 

the coefficients would have been related as shown in (4. 25) and the following 
expression for the field would have been obtained: 


E h ®> 


-itdju.C 

lee 


47 r 


ikjR 

R 



(4. 37) 


Equations (4. 36) and (4. 37) may be compared to the expression for the far field 
radiated from an ideal Huygens' source in free space (4. 11), repeated here: 


e h Cr> = 


-iu/i C 
lee 


47T 


ik R 


cos 0) cos ^0 + (l - cos 0)sin^ ^ 


(4. 38) 


It is seen from (4. 36), derived using = rj^C , that even a very small sphere 

has an effect on the radiation pattern. In particular, the contribution from the 

electric dipole source has an additional factor of 3/(e^+2) as compared to the 

free-space case. The magnetic source contribution, however, has remained 

unchanged. Note that when e = 1, equation (4.36) reduces to (4. 38), the 

r 

free-space case, as it should. 

The small dielectric sphere apparently behaves like a short electric 
dipole with a current moment 


C 

e 


_3 

e +2 
r 




1 - e 
r 

e +2 


r 


(4. 39) 
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which is independent of the sphere size. This factor also appears in the results 
found by Stratton (1941, p. 206) for a dielectric sphere within a uniform static 
electric field. 

Let us now look at the deviation of the small dielectric sphere patterns, 
using C *n, C or C = rj n C for the electric and magnetic current relations, 
from the pattern radiated by a Huygens' source in free space. 

Let 


O'. 


e +2 
r 


and 


O' 


!£ 
e +2 
r 


The first terms of (4. 36) and (4. 37) then contain the respective factors 

1 - <*jCos 0 

(1- accost?) 

and similarly in the second terms we find the factors 



- cos 0 


(a 2 - cos« . 


For the Huygens' source in free space a = a - 1 and a (1 - cos 0) pattern 

1 £t 

remains. 

Table 4-1 shows the values of a and for various ^values of dielectric 
constants. The a ^ remains a good deal closer to unity than o^, and hence 
using rjg for the current relations in constructing the Huygens' source on the 
surface of a small dielectric sphere will yield only a small error. In fact, 
differs from unity by only 6 percent at its maximum value corresponding to a 
dielectric constant of 2. 
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TABLE 4-1 


Approximate Huygens' Source Factors 
For Small Spheres 


e 

r 

*1 

a 2 

1.0 

1.0 

1.0 

1.5 

0.86 

mem 

2.0 

0.75 

mm 

2.5 

0.67 

1.05 

3.0 

0.60 

1.04 

3.5 

0. 55 

1.02 

4.0 

0. 50 

1.00 

4.5 

0.46 

0. 98 

5.0 

0.43 

0. 96 




Chapter V 


PROPERTIES OF A HUYGENS’ SOURCE 
IN THE PRESENCE OF A DIELECTRIC SPHERE 


5. 1 Introduction 

In this chapter, the orthogonality relationships between the n and in 
functions are shown and then an expression for the radiated power is derived. 
Equations for the radiation resistance and directivity are formulated and com- 
pared to the free-space case for small spheres. Finally, an expression is 
developed for the power dissipated due to the nonzero conductivity of the dielec- 
tric material. 

5. 2 Radiated Power 

In the far zone we found (4. 30, 4. 31) that the electric field radiated by 
a Huygens' source in the presence of a dielectric sphere may be represented by 


e h ( R) 


-up.c 

lee 


ik,R oo 


47T 


R 


— - + 0- n ~ i -> 


•1 


n=l 


n , — 

' m 

n oln n eln 


(5. 1) 


where 


a = -i(A +D ) , 

n n n 

g =B (1) + 0 (1) 
n n n 


or 


a = -i[A +D ) , 
n n n 

g =b (2) + c (2) 

n n n 


for the source located in Regions I or II respectively. The vector wave functions 
in (5. 1) are found from (4. 28) and (4. 29) to be 


m 


oln 


P*(cos 0) 
n 

sin0 



BP'*’ (cos 6) 
n 

ae 
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3pl(cos 0) 
n 

30 


COS I 


(cos 0) 
n 

sin 0 



(5. 2) 


An expression for the total radiated power W may be derived by 
integrating the total power flow across an infinitely large sphere, centered at the 
origin, using the relation 


* 27 r nir 


W = 


2 n 


vly 2 

e • e'r dn 


(5.3) 


0 


>0 oo 


where dO= sin 0d0d^ , and * denotes the complex conjugate. 

Let us now investigate the various scalar products which will be formed 

in (5. 3). From (5, 2), it can be seen that each term of the product m , -n , . 

' oln el* 

will contain a factor of the form 

p!(cos0) 3pl(cos0) 

n l 

sin0 30 


Using the differential relations for the associated Legendre functions (Stratton, 
1941, p. 402), it is found 


P 1 (cos 0) 3pl(cos 0) 
n l 

sin0 30 
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sin0 


P^ (cos 0) - 


P (cos 0) 
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sin0 


P f (cos 0) I 


It now follows that 


pl(cos0) 3pl(cos0) 
n i 


sin0 


30 


sin0d0 = 0 


(5.4) 


since 



Thus 


m . *n . 6 
oln elf 


dfi = 0 . 


Also, on performing the integration. 
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m ‘ m = \\ n • n . dfi = 27m (n + 1) g 

oln olf \\ eln ell 2n + l ni 


(5.5) 


where 


fO, nf 

6 nf " { 

LI, n = 


and we have used the formula (Stratton, 1941, p. 417) 

/dP m dP™ p m P m \ 

+ 2 P n P i \ . , in 2n(n-H)(n+m)! , 

, n \ d0 d0 . 2 J 1 (2n + lXn-m)! ni 

;0 \ sin 0 / 

(5. 6) 

Upon inserting (5. 1) into (5. 3) and evaluating the integrals with the aid of (5. 5), 
we obtain the desired relation for the radiated power 


n=l 


oo 


(5.7) 


where r) Q = 1207 rft. 


5. 3 Radiation Resistance ( 

Suppose that instead of having an infinitesimal point source, the electric 
current element were of a short length i with a constant magnitude 1^. Then 
C g = I Q i, and the radiation resistance for this system is easily obtained. 

The radiation resistance of an antenna is defined as 
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That is, assuming no losses, a current element with an rms value of 1^ flowing 
on an antenna with a radiation resistance R will radiate an average power W. 

The far-zone electric field of the crossed dipole source in the absence 
of the dielectric sphere is given by (4. 10). Using this equation together with 
(5. 3) and (5. 8), the radiation resistance is given by 

R = 20(k 1 i) 2 fl + ^M . (5.9) 

It is apparent from this equation that for the Huygens' source in free space, i.e. 

2 

e =1, the radiation resistance obtains its maximum value of 40 (ki) which is 
r 

twice that of an elementary electric dipole source. 

To determine the radiation resistance of the source in the presence of 
the dielectric sphere (5. 7) is substituted into (5. 8). Hence 

oo 

R = 15(k,i) 2 / . (2n + l)(a,a*+/3j3*) . (5.10) 

1 * — rf n n n n 

n=l 

Let us now place the source on the surface of the sphere, allow the 
radius to become small, and see how the resulting radiation resistance com- 
pares to the free-space case. Using the relations (4. 35) to evaluate a and )3 

we find for the case when C = n,C 

m 1 e 


R = 20{k-0‘ 


t-m 


(5. 11) 


Had we let C = n 0 C in relating the electric and magnetic current 
m 2 e 

moments, instead of (5. 11) we would have obtained 


R = 20 (ki)' 


[r*W] 


(5. 12) 


In order to compare the radiation resistances for the various sources, (5. 9), 

2 

(5. 11) and (5. 12) are divided by (ki) and the results denoted by R , R and R , 

U -L Z 
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respectively. The results are shown in Table 5-1 for various values of e . 


TABLE 5-1 

Normalized Radiation Resistance 
for Sources on Small Spheres 


m 

R o 

R i 

R 2 

i 

40.0 

40.0 

40.0 

2 

30.0 

31.2 

21.2 

3 

26.6 

27.2 

13. 8 

4 

25.0 

25.0 

10.0 

5 

24.0 

23.6 

6.6 


It can be seen that the normalized radiation resistance R^ of the source 
(C^ = 17 ) on sur ^ ace °f a sn ?all dielectric sphere closely corresponds 
to the free-space case, R^. It should be noted, however, that R^ does not cor- 
respond to the Huygens' source, i.e., l-cos0 pattern, unless = 1 in (5.9). 

The normalized radiation resistance R„ of the source (C = n n C ) is smaller 

z m Z e 

than R^ because in the former the magnetic current moment is reduced in order 
to compensate for the reduction of the electric current moment by the dielectric 
sphere. In the next chapter, it will be seen that the radiation resistance for 
larger spheres tends to oscillate about the free-space value. 

5.4. Directivity 


The directivity in the forward direction relative to an isotropic source 
is defined by 


_ 47r(radiation intensity in 0 = 180° direction) 
total power radiated 


(5. 13) 


since for the sphere concentric with the origin, and the source located on the +z 
axis, we are interested in the relative field intensity in the -z (0 = 180°) direction. 
From (5. 2), it is evident that 
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and 
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n+1 n(n + l) 


6=n 
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90 
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(5. 14) 


Using this result in (5. 1), an expression for the radiated power intensity 
in the -z direction is obtained. Thus, 



— — Jr 
E -E 


6=ir 


2 2 
15 k C 
1 e 

47T 
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2n + l 
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(i) V +j3 ) 
n n 


(5. 15) 


where rj \2tiir£l. 

Combining (5. 15), (5.7) and (5. 13), the desired directivity formula 
results : 


D = 


e 

n=l 

oo 

s 


(ilV+ft) 

2 n n 


2n + l t *^0 
— — (a a +f3 p ) 
2 n n n n 


(5. 16) 


Again, it is noted that A = B = 0 or C = D = Oin the expressions for a 

n n n n n 

and in Section 5. 2, in order to recover the result due to the magnetic or 
electric source only. 

It is found from (5. 16) for the small sphere case that 
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and 


D„ 




2 

) +1 


for C 

m 


"2 C e 



The results are shown for various values of e in Table 5-2 . 

r 


TABLE 5-2 

Directivity of Sources on Small Spheres 


e 

r 

V 

°2 

1 

3.000 

3.000 

2 

2. 940 

2. 997 

3 

2. 824 

2. 999 

4 

2.700 

3.000 

5 

2. 586 

2. 999 


(5. 17) 


(5. 18) 


From (4. 11) and (5. 13) it can be shown that the directivity of a Huygens' 
source is 3. From Table 5-2 it is seen that the directivity D 2 remains 
approximately this value. 

In Chapter 6, it will be of interest to compare the directivity of the 
pattern radiated by the dielectric sphere with that of a uniformly illuminated 
aperture of the same diameter. The normalized field pattern of a large uniform 
aperture as shown in Fig. 5-1 can be found by the application of Huygens' princi- 
ple and is given by Kraus (1950): 


E(0) = 


2X 

7rd 


J i(f sini j) 

sinS 


(5. 19) 
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uniform 

plane wave 


T 

d 

_L 




infinite sheet 


FIG. 5-1: UNIFORMLY ILLUMINATED APERTURE. 


For this case, the expression 


D = 


E(0)P 


^|E(e)| 2 df2 

is used for the directivity, and hence 


D = 


2.2 

TT d 


2X 
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7T/ 2 [j t (x 


sin0 


d0 


5.5 Dielectric Loss 


For a lossy dielectric e 2 * s complex and may be represented by 
e 2 = e’+ie” , 


where e' and e" are the real and imaginary parts of the permittivity 


(5. 20) 


(5.21) 



52 


respectively. The conductivity <j is related to e" by 


a = toe" , 


and the loss tangent is defined as 

tan6=e ,, /e , . 

Equation (5. 21), therefore, may be cast in the form 


Cg = e’U+itanS) . 


The propagation constant k will also be complex, being equal to 


k = k'^l+itan 6 . 
6 


Thus, in the previously derived equations, k^ and also are simply 

replaced by their complex values to account for the dielectric loss. 

To find the power dissipated within the sphere, the source is placed 

in Region I close to the surface of the sphere. With C = r) n C , it is found 

m i e 

from (3. 6) and (3. 18) that the electric field within the sphere is given by 


h® = {‘' M 0 i tl (k 2 > + 5N elr , (k 2 >} • 

n-i J 


“oiA* -J.W 


P (cos 0) 


. 9P (cos 0) 

rZ n 


>in^ % 


N . (kj = (p)p^( C os 0) cos ^R 

ein 2 p n n 

, rap 1 (cos 0) 


P (cos 0) 
n 

sin0 


(5. 28) 


P=k 2 R 
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( 1 ) 


m ic [p h (p )"]’ 

(1) v, u — S- L. 3 a 3J 


v = c h (p ) - 
n n 3 




,(D 


fp„h (p fl« d' 
; =d «) L 3n . 3j +1 n h 


n P, 


r “ \ 

F “ 


, to ,(D . . „ ... 

and c , d are given by (2. 14). 
n n 


The power dissipated within a region of conductivity a is given by the 


formula 


P D = 2'' ' aE ‘ E dV * 


(5.29) 


For a homogeneous sphere of radius a (5. 29) becomes 


p27rp>7r na 
P D “ \ 

do do UO 


E (R) ‘ E (R) R 2 sin 9 dR d6 dj£ . 


(5. 30) 


Because of the orthogonality of the associated Legendre functions and 
their derivatives, shown in Section 5. 2, 

M • N , „ sin 6 d0 s 0 . 
oln eli 

Jo 

Also, it is easy to show 
'27r rw 


— * 


M oln W ' M ol< (p,dn= 2 ™ 2 


JO 


and 


\2-n f\n 


Jo JO 


N„ 1 _(p) • N„ n 0 (p) df2 = — ir ~ + ^ [n(n + l)j n (p)j J(p) 


eln lx eir 


pp (2n+l) 


[pj^' , (5.31) 
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where we have used (5. 6) and the formula (Stratton, 1941, p. 403) 


P m (cos0)p“(cos0)sined0 = 7 — 7 6 . 
n i 2n + 1 (n - m) ! a? 


With the aid of (5. 31) the 6 and 0 integrations are performed and an expres- 
sion for the dissipated power obtained. Thus, 


2 2 2 2 oo 


+ n(n + l)i (p)j*(p) + 


n ' n 


[pj n <P)]' [R 2 dR 


15 ot) ^ , pa 

P D = ~ [ r j~ e E <2a + l)\ ||^ n (p)| 2 + n(n + l)|5i n (p)| 2 

+|l[pj n (p)]'| 2 \dR . (5.32 

Now that expressions for the dissipated power P^ and the radiated 
power W (5.7) are available, the antenna efficiency may be computed from the 
relation 


Power radiated W 

Total input power W + P^ 


(5. 33) 



Chapter VI 

NUMERICAL RESULTS 


6. 1 Introduction 

Some of the more significant equations derived in Chapters IV and V 
were programmed and evaluated using the University, of Michigan EBM System 
360 Model 67 dual-processor computer. The data was then plotted on the 
CALCOMP 780/763 digital plotter. 

In this chapter, the numerical results using a Huygens' source are first 
compared with measured data for plexiglas spheres placed over the ends of 
waveguides. Curves are then presented showing the variation of directivity with 
sphere diameter for various dielectric constants. Subsequently, the resonant 
modes of the dielectric sphere are studied together with the effect of dielectric 
loss. Finally, curves are presented which show that directive patterns may be 
produced with the source located inside the sphere. 

6.2 Comparison with Measured Data for the Plexiglas Sphere 

The far-zone electric fields radiated from a Huygens' source on the 
surface of a dielectric sphere (4. 31) were calculated numerically for the two 
principle planes using the computer program listed in Appendix B. Since the 
series for the far field converges in a manner similar to the spherical 
Bessel functions it was found that N = kD +7 terms provided sufficient accur- 
acy. 

The normalized power patterns for various diameter spheres, computed 
in one degree increments, are presented in Figures 6-1 to 6-4. Here the E and 
H planes correspond to the = 0 and / = 7r/2 planes, respectively. In each case, 
a dielectric constant of 2. 57 with a loss tangent of . 0065 was used., which is 
typical of plexiglas in the frequency range of interest (von Hippel, 1954). It 
should be noted that with the source oriented as shown in Fig. 4-3, the maximum 
radiation occurs in the -z direction, and hence the top of the computed patterns 
correspond to Q = 180 degrees. 
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The curves on the left of each figure were measured by Croswell 
and Chatterjee in the anechoic chamber located at the NASA Langley Research 
Center. These measurements were made using plexiglas spheres centered 
against the flange of an open-ended waveguide. The parameters are listed in 
Table 6-1. It should be noted that no attempt was made to impedance match the 
waveguide or, in the case of the rectangular guide, to physically fit the aper- 
ture to the curved surface of the sphere. 


TABLE 6-1 

Measurement Parameters 


Sphere Diameter (inches) 

3 

5 

5 5 

Frequency (GH ) 

5.0 

5.0 

8.0 10.0 

Sphere Diameter (X Q ) 

1.27 2.12 

3.39 4.23 

Waveguide Type 

C-band 

Circular 

X-band 

Rectangular 


From these patterns. Figs. 6-1 to 6-4, it can be seen that the agree- 
ment is quite good, particularly for the larger diameter spheres. The theo- 
retical patterns do, however, show sharper nulls in the main beam and stronger 
back lobes than the measured patterns. These differences may be due to scat- 
tering from the waveguide feed structure used in the measurements. 

Figure 6-5 shows the directivity of the Huygens' source on the surface 
of the sphere computed from (5.16). Again, a dielectric constant of 2.57 and 
loss tangent of . 0065 have been used. For comparison, the directivity of a 
uniformly illuminated circular aperture (5. 20) is also shown. The abscissa 
is measured in aperture diameter or sphere diameter in free-space wavelengths. 
The small circles and crosses in the figure represent directivities derived 
from measured patterns for waveguide-excited plexiglas spheres using the 
empirical formula 

K 


D = 
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EXPERIMENTAL (E-PLANE) THEORETICAL (E-PLANE) 


Huygens’ Source 


€ ■ 2.57 
r 

tan 6= . 0065 



EXPERIMENTAL (H-PLANE) • THEORETICAL (H- PLANE) 


FIG. 6-4: RADIATION PATTERN FOR PLEXIGLAS SPHERE (D = 4.23A q ) 
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FIG. 6-5: DIRECTIVITY OF PLEXIGLAS SPHERE (e = 2.57, 
tan6 = .0065). r 


where K is a constant and 9 and are the half- power beamwidths in the 
two principle planes measured in degrees. By comparing the computed directivity 
with the beamwidths of the Huygens' source patterns, it was found that K = 26, 500 
yielded an error less than 0.5dB for sphere diameters in the range of interest. 
Since the measured and computed patterns have the same general shape, this 
value of K was used in the above formula. 

It can be seen from Fig. 6-5 that the experimental and theoretical 
directivities agree quite favorably with a maximum deviation of 1 . 7dB occurring 
at D = . 85 It is also noted that for the measured data available, the direc- 
tivity is generally greater than that for a uniformly illuminated aperture. 
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In Table 6-2, the radiation patterns presented in Fig. 6-3 for the 
dielectric sphere are compared to that observed for a conical horn of about the 
same diameter (3. 4X) as given by Jasik (1961). The data is for a horn opti- 
mized in the sense that the flare angle has been adjusted to yield a maximum 
phase deviation of 3/8X in the aperture wavefront. 


TABLE 6-2 

Comparison with Optimum Horn ' 


Parameter 

Optimized Horn 
D = 3.4 X q 

Plexiglas Sphere 
Measured 

D = 3. 39X q 
T heoretical 

E -plane, 3dBbeamwidth 
(degrees) 

16 

12. 5 

14 

H-plane, 3 dB beamwidth 
(degrees) 

19 

14 

14 

First sidelobe level, 
E -Plane (dB) 

- 8 . 5 

-13. 5 

-16 

First sidelobe angle from 
axis. E-plane (degrees) 

.22 

20 

22 

Directivity (dB) 

17.7 

23.7 

21.0 


From this table, one can see that both the measured and theoretical 
patterns for the plexiglas sphere have lower sidelobes and higher directivity 
than the optimized conical horn. 

Figure 6-6 shows the radiation resistance computed from (5. 10) toget- 
her with power reflection coefficient curves measured for 2 and 3-inch plexiglas 
spheres. 

Here it is noted that both the radiation resistance and the measured 
curves have strong peaks occurring with a periodicity of about D = 0.25X^ or 
(k 2 a) =1.26. 

The antenna efficiency was computed using (5.33) and is plotted in Fig. 
6-7. For this computation, the source was placed O.OIX^ from the sphere 
surface to avoid the problem of infinite loss associated with a point source in 
contact with a lossy medium (Tai, 1947). 
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As one might expect, the antenna efficiency generally decreases with 
increasing sphere diameter. It is also noted (Fig. 6-7) that for the larger 
sphere diameters the efficiency curve shows a marked decrease at the same 
points the radiation resistance (Fig. 6-6) shows sharp peaks. 

6. 3 Directivity 

Since the directivity is an important parameter in antenna design, 
the value of dielectric constant which maximizes this parameter for a given 
diameter is of interest. 

The directivities for the Huygens' source on the surface of the sphere 
(5. 16) are plotted as a function of sphere diameter for various lossless dielec- 
trics in Fig. 6-8. Since the curves were plotted at only a finite number of 
points (0.01 diameter interval), it is recognized that some of the fine structure 
may not be precise. Nevertheless, it is apparent from these figures that as the 
diameter increases, the directivity becomes frequency-sensitive, with the effect 
being more pronounced for the higher dielectric spheres. The sharp dips in 
these curves are due to a resonance effect, discussed in detail in Section 6.4, 
caused by the dielectric sphere in which much of the energy is radiated into 
the sidelobes. 

The resonance effect was not noted in the relatively smooth curve of 
Fig. 6-5 because the dielectric loss was included in the calculation. By com- 
paring the directivity curves of Fig. 6-9 which were calculated using a loss 
tangent of .0065 with the lossless curves of Fig. 6-8, it is noted that the loss 
tends to attenuate the resonance with the effect being most notable at the larger 
sphere diameters. The dielectric loss, however, had little effect on the peaks 
of the curves. 

Figure 6-9f was constructed by drawing smooth curves through the 
maximum points on the directivity curves presented in Figures 6-5 and 6-9. 

Here it is seen that a dielectric constant of about 3. 0 tends to optimize the direc- 
tivity, particularly for the larger diameter spheres. 



DIRECTIVITY (DB) 

5.0 10.0 15.0 20.0 25.0 


.5 1.0 l.S 2.0 2*5 3.0 3.5 >4.0 4.5 S.O 

DIAMETER IN WAVELENGTHS 



FIG. 6-8: DIRECTIVITY OF LOSSLESS SPHERES. 
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FIG. 6-8: DIRECTIVITY OF LOSSLESS SPHERES. 



'•0 -s 1.0 1.5 2.0 2^5 3lo 3^5 4.0 4.S 5.0 

diameter in wavelengths 

FIG. 6-8: DIRECTIVITY OF LOSSLESS SPHERES. 





'.0 .5 1.0 1.5 2.0 2.S 3.0 3.5 4.0 4.5 5.0 

01RHE7ER IN WAVELENGTHS 

FIG. 6-9: DIRECTIVITY OF SPHERES WITH TAN 6 = .0065. 


DIRECTIVITY (OB) DIRECTIVITY CDB) 

1 0.0 15.0 20.0 25.0 .0 5.0 10.0 15.0 20.0 25.0 
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(d) e = 4; 0 
r 


•-0 .S 1.0 l.S 2.0 2.5 ■ 3.0 3.5 4.0 4.S 5.0 

0IRME7ER IN WAVELENGTHS 


FIG. 6-9: DIRECTIVITY OF SPHERES WITH TAN6 = .0065 
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6. 4 Resonance and the Effect of Dielectric Loss 

The resonance phenomenon noted in the directivity curves of Section 
6.3 is now investigated. Figure 6-10 shows a detailed section computed using 
a diameter interval of . 003 X. A pattern occurring in this curve is noted. For 
example, resonances occur at B, D, and then the pattern repeats starting at E 
with the nulls becoming stronger with increasing diameter. 

Figure 6-11 shows the antenna patterns corresponding to points B, C 
and D of Fig. 6-10. The patterns of Fig. 6-lla (H-plane) and Fig. 6-llc 
(E-plane) each have ten complete lobes indicating that the resonance is due to a 
large tenth mode. It has been found that points A and F are caused by 
resonances in the 8th and 13th modes respectively. Hence, with increasing 
diameter, the resonance appears in the H-plane, then the E-plane, and 
again in the H-plane at the next higher mode. 

Looking back at (5. 1), it is seen that the far-zone field in the two planes 
may be expressed by 


E(R) = 


_ ik.R co 

^l C e e 1 V 2" + l , •>» 

— Z-f <-*> 

n=l 


Air 




P (cos 6) 


V. 


r p (< 
U-2L 

I n si 

r ap 1 

\a— a 


sin0 
1 


+ 0 . 


n(n+ 1) 


9P (cos 0) 
n 


n 


+/3 


90 


P (cos 0)* 
n 

sin0 


A 

e , 


t. 


E- plane 


H-plane 


( 6 . 1 ) 


It has been found that when there is an ith mode resonance in the E-plane or 
H-plane the magnitude of j3^ or a., respectively, is considerably larger than 
the other coefficients. It therefore can be seen from (6. 1) that the radiation 
patterns will exhibit a strong resemblance to 


C 


ap.Vos 0)\ 2 


90 


- 


in the plane of resonance , 



. X 1 1 1 1 

1.5 2.0 2.5 3.0 3.5 

DIRMETER IN HRVELENGTHS 

FIG. 6-10: DIRECTIVITY OF HUYGENS’ SOURCE ON SPHERE 

SURFACE (e = 3.2, tan 6 = 0.0) 
r 
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E-PLANE (c) 0*2. 54 A H- PLANE 

o 


FIG. 6-11: RADIATION PATTERNS FOR VARIOUS SPHERE 

DIAMETERS (e =3.2, tan 6 =0.0). 
r 


Pj ( cos 9 )V* 
sin 9 J 


in the alternate plane. 


This resonance phenomenon occurs when the denominators of (3. or cr , 

which contain the respective factors e r |p^h n ( )j ' j^Pg) ~ [ P 2^n^2^]' ^P^ 

or [f^n P 2^] ' ^n ^1^ become small. It has been shown by 

Stratton (1941, Sec. 9.25), however, that these factors never reach zero. 

The total radiated power may be expressed as a summation over the 


power radiated in each mode, i. e. , 
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CD 



P 

as shown by (5.7). The quantity — • 100 percent represents, therefore, the per- 

w 

cent of the total power contributed by each mode. These power contributions are 
plotted for the modes with significant power for points B and C in Fig. 6- 12a 
and F and G in Fig. 6- 12b. Although the modes are discrete they have been 
connected for identification purposes. 

From these figures it is noted that when the directivity is high, the power 
is distributed among the modes. At the points of resonance, however, most of 
the power is concentrated in a single mode; about 89 percent is in the 10th mode 
at B, and 98 percent is in the 13th mode at F. 

Although it is not apparent from Fig. 6-12 because of the normalization, 
the total radiated power at resonance is considerable higher than at nonresonance, 
and hence the radiation resistance (see Fig. 6-13) shows strong peaks at the reso- 
nant points. This curve was computed using (5. 10) in .001 diameter intervals. 

It was noted from comparing the directivity curves for lossy spheres 
with the corresponding lossless curves that the attenuation had little effect on 
the points of high directivity. The dielectric loss did, however, significantly 
reduce the magnitude of the resonances with the effect being more pronounced 
for the larger diameter spheres. The reason for this may be explained with the 
aid of Fig. 6-14. 

We define modal efficiency by 

P (radiated) 

Efficiency = p (radiated) +P (loss) ' 100 1,61:06,11 ' 
n n 

where P (radiated) and P (loss) are the nth terms of (5. 7) and (5. 32) respec- 
n n 

tively. In Fig. 6-14 this efficiency, using a loss tangent of .01, is plotted as 
a function of mode number for various resonant sphere diameters. Generally, 
as the sphere becomes larger, it is seen that the efficiency of the higher modes 
increases while that of the lower modes decreases slightly. 



Relative Power (percent) 
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FIG. 6-12a: RELATIVE POWER DISTRIBUTION AMONG MODES. 



Relative Power (percent) 
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FIG. 6- 12b: RELATIVE POWER DISTRIBUTION AMONG MODES 
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FIG. 6-13: RADIATION RESISTANCE FOR HUYGENS' SOURCE 
ON SURFACE OF DIELECTRIC SPHERE (e =3.2, 
tan 6 = 0.0). r 
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RESONANT FREQUENCIES. 


Efficiency (percent) 
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FIG. 6-15: MODAL EFFICIENCY FOR VARIOUS LOSS TANGENTS 
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Table 6-3 shows that each new resonant mode has an efficiency less 
than the previous one in the same plane. This, together with the fact that for a 

TABLE 6-3 


Modal Efficiency for Resonant Spheres 
(loss tangent = . 01; = 3. 2) 


Diameter (A^) 

Mode No. 
Number 

Resonance 

Plane 

Efficiency (o/ ) 

'O 

2.000 


8 

H 

48. 55 

2.407 


10 

H 

25.22 

3.010 


13 

H 

6.08 

2. 110 


8 

E 

61. 61 

2.533 


10 

E 

36.78 

3. 145 


13 

E 

9.48 


given diameter the resonant mode efficiency is less than that of the lower modes, 
accounts for the smoothing of the directivity curves with dielectric loss . These 
results also explain the reduction of the overall antenna efficiency, noted in 
Fig. 6-7, in the resonance region for the larger sphere diameters. The effect 
of other loss tangents on the modal efficiency is shown in Fig. 6-15. 

The patterns of Figs. 6-16 and 6-17 were computed in 2-degree intervals 
using (4. 13). Because of the high loss used in Fig. 6-16d, the internal reflec- 
tions are greatly attenuated, and hence there is practically no radiation in the 
back direction. Finally, for a large sphere with constitutive parameters typical 
of plexiglas, Fig. 6-17c, the resonance is almost completely attenuated and 
therefore may not be noted in actual measurements. 

6. 5 Source Displaced from Sphere Surface 

To examine the effect of moving the source away from the sphere, the 
curves of Fig. 6-18 were calculated using (4. 30) in two-degree intervals. The 
dielectric constant was kept at 2. 57 with a loss tangent of .0065. Since b is the 
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E-PLANE (c) tan <5 = .0065 H-PLANE 


FIG. 6-17: RADIATION PATTERNS FOR SPHERES WITH 

VARIOUS LOSS TANGENTS (D = 4.21 X , 

e = 2.57). ° 

r 

distance from the origin to the source and a is the sphere radius, the ratios 
of b/a shown in Figs. 6-18a, b, c and d correspond to the source spaced 0.0, 
0.25, 0. 5 and 1.0 free-space wavelengths from the sphere surface. From 
these patterns it is seen that the main beam first narrows slightly and then be- 
comes broader with increasing displacement. The sidelobe level also becomes 
considerably higher as the source is moved away from the sphere. 

Without varying the dielectric constant, the curves shown in Fig. 6-19 
were computed using (4. 31). These patterns correspond to the source spaced 
0. 1, 0. 5, 0.75 and 1. 5X^ inside the sphere surface. The theoretical patterns 
shown in Fig. 6-3 were also calculated using (4. 31) for the source on the sphere 
surface, and may be compared to Fig. 6-18a which was computed from (4. 30) 
using the same parameters. 



87 


From Fig. 6-19, it is seen that the antenna pattern continuously de- 
grades as the source is moved into the sphere, and that the sphere shows very 
little focussing effect as the source approaches the center. 

Directional patterns can, however, be obtained with the source located 
inside the sphere by appropriately adjusting the dielectric constant as shown in 
Figs. 6-20 and 6-21. The directivity of these patterns is about the same as 
that with the source located on the surface of a sphere with a lower dielectric 
constant. 

Figure 6-22 shows the optimum ratio of source to sphere radius for 
various diameter spheres. These points were obtained by fixing the dielectric 
constant at 4, 6, 9, and 12 and varying the source position until the maximum 
directivity was found. The additional points corresponding to the source on the 
sphere surface were obtained by adjusting the dielectric constant. In each case 
a loss tangent of . 0065 was used. For comparison, the upper curve shows the 
location of the geometrical optics paraxial focus (8. 5) . 

It is seen from Fig. 6-22 that for a dielectric constant of 9 the source is 
located about half way between the center and back surface of the sphere. It is 
also noted that for a fixed source position a higher dielectric constant is neces- 
sary to obtain optimum directivity for the smaller spheres. In all cases, how- 
ever, the best results were obtained for a dielectric constant less than that 
required for the paraxial focus. 







E -PLANE 


(d) *-1.59 


H-PLANE 


FIG. 6-18: RADIATION PATTERNS FOR SOURCE LOCATED 
OUTSIDE SPHERE (D = 3. 39 X , € = 2.57, 
tan 6 =.0065). ° r 









E-PLANE 



FIG. 6-20: RADIATION PATTERNS FOR SOURCE LOCATED 
INSIDE SPHERE (D = 3.39X , tan 6 = .0065). 





FIG. 6-22: SOURCE LOCATION AS A FUNCTION OF DIELECTRIC 
CONSTANT FOR OPTIMUM DIRECTIVITY. 



Chapter VII 
EXTENDED SOURCES 

7. 1 Introduction 

The previous chapters were concerned with the effect of a dielectric 
sphere on the radiation from infinitesimal point sources. In this chapter the 
excitation is expanded into line sources concentric with the sphere. 

Expressions for the far field radiated by curved electric dipoles are 
derived first, and then the results for a superposition of Huygens' sources in 
the presence of a dielectric sphere are presented. Finally, it will be shown 
that the patterns from these sources do not differ greatly from those of point 
sources. 

7.2 Electric Field Due to a Curved Electric Dipole 

Consider a dielectric sphere with a filament of current located in the 

x - z plane at R = b as shown in Fig. 7-1. It should be noted that in the ^ = 0 

z 



FIG. 7-1: CURVED FILAMENT OF CURRENT IN PRESENCE OF 
DIELECTRIC SPHERE (x-z PLANE). 
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/ A A 

and p = 7T half-planes the current is defined to flow in the +0 and -0 directions 

respectively. This is because the 0 coordinate variable is always measured in 

a downward direction from the +z-axis . 

For the antenna located as shown in Fig. 7-1 we let 


•»»'> = I(6,) [}<!») -#*'-»»] S ’ 


<7. 1) 


where 1(0*) is the current distribution on the antenna in the j4 = 0 half-plane 
and is only a function of the coordinate variable 0*. 

To find the electric field in Region I, (7.1) is substituted into (2. 1). Thus, 


E(R) = iwju. 


' d 0 f' 1r n oo 


JO ^0 dQ 


g“(r|r')*0' i( 0') 6 — * --- b) 


[s$') - 7rj] R* 2 dR'dji'd©' . 
Integration over the R’ and variables yields 

’o 


(7.2) 


E(R) = iuju^b 


KOOG^CrIr') 


R ,= b 


N0, 


J =n _ _ 

I(0’)G 3 (R R') • 0'd0* 


R’=b 

j) x -TT 


(7.3) 


Following the same procedure for equations (7. 1) and (2. 2), a similar 
result is obtained for the electric field in Region n. Thus, 


\0 


E(R ) = iu/u^b 


I(0')G 21 (R R'J-D'd©' 


R'=b 

^ ,= 0 


'0 


I(0')Go 1 (R | R ') • 0 f d0' 


(7.4) 


R'=b 
j> x =i T 
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The expression G^(R jR') ' 0 must now be examined. From equations 
(2. 8) and (2. 9) it is seen that at R‘ — b 


M* (k) • 0’ = N' (k) • 0 ' = 0 
emn omn 


U' = 0 or 


__ ^ P m (cos0') 

M' (k) • 0' = mj (p ) ; "\ 

omn n 3 smfl' 


- [ p 3 i o (p 3^' 

N' «.«• = ^ — y 

emn p„ 90* 


>*' =o 


(7.5) 


M’ (k) • 0 f 
omn 


' = 


= (-l) m M' (k) • 0' 

omn . 

IT p'=0 


N' (k) • 0’ 
emn 


= (-l) m N' Oc)-0’l 

rf'-jr emn 1^=0 


(7.6) 


where p denotes kb. 

O 1 

Substitution of (2. 16) through (2. 18) into (7. 3) and (7. 4) yields the 
following expression for the electric field in the various regions: 


e(r )= Y fa-sat ; . 

JLW 7t n(n + 1) (n + m)! 

n=l m=l 

Xm ( 1) (kJA (1) I (0 )+N (1) (kjB (1) I* (0jl R > b > a 
(_ omn 1 n mn 0 emn 1 n mn 0 J 

(7.7) 
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(7. 10) 


and 2^ is used to denote that the summation is over the odd integers only. The 
scattering coefficients are given in Section 2. 2. 

The far-zone field of (7. 7) is easily found with the aid of (4. 27). Hence 
in this region 


E(R) = 


-u/^b 


7T 


ikjR 


R 


cp n 

£2 

n=l m=l 


2n + 1 (n - m) ! 


n(n + l) (n+m)! 


(-i) 


n 


f-im A (1) I (©J+n B (1) I' (©_)!-, (7.11) 

omn n mn 0 emn n mn 0 I 

where and are given by (4.20). 

For computational purposes, it is often desirable not to have the 
associated Legendre functions within the integrals of (7, 10), It will now be 
shown that by locating the source in the x - y plane, the integration will be over 
the variable^ and hence the integrand will contain only trigonometric functions. 
Consider a curved filament of electric current oriented as shown in 


Fig. 7-2. In this case, let 
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FIG. 7-2; CURVED FILAMENT OF CURRENT IN PRESENCE OF 
DIELECTRIC SPHERE (x-y plane). 
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jCr') = i(^') 5(0'- f)? 1 


and assume that the current distribution l(f l ) is an even function of the coor- 
dinate variable (' . 

To solve for the fields in Region I, (7. 12) is substituted into (2. 1) and 
on performing the R' and 0' integration one finds 


E(R) = iu/Ltb\ Gg 1 (R|R , )-£l({* , )d^ , | 


R'=b 
0 , = 7t/2 


+iUA*b \ G^ 1 (R|R')‘?I(^ , )d^| 


R'=b 
0 , = 7r/2 


Alternately, since I(jJ') is an even function, this result may be expressed 
in the form 


E(R) = 2iujub \ Seven^ 1 l R ^ * ?* \ 


R*=b 
0 ,= 7t/2 


where is an even function of j>\ Taking the even parts of the spher- 

ical wave functions (2. 8) and (2. 9) results in 

A A 

M t .ft = N I .ft = o 
omn r emn 

311(1 _ A 9P m (0) 

M* • 0' = -j (p ) — ^ — cos m^' , 

emn , n 3 90 

R ,= b 

j>' = ir / 2 


- S) 
N* -A 
omn 


["aWET m, . 

mP (O)cosmp' 


R'=b 
^’ = 7r/2 
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It is known (Smythe, 1950, p. 153) that 

n-m 


p m (o)= s 

n 


C , 2 1 < 3 • 5. . . (n+m - 1) 

2 • 4 • 6 . . . (n - m) 




(- 1 ) 
0 , 


and also 


3p“<0) m+1 

~t~ ~ - P n <°> ■ 


n + m even 


n + m odd 


Using these relations in (7. 14) the result is obtained 


where 




n=l m-0 


n(n+ 1) *m^0^ 


M (1) (k 1 )A (1) P m+1 (0) , 
emn Inn 


l N (1) (k )B (1) mP m (0) 
v omn In n 


n + m odd 


n + m even 


R>b 


(7. 16) 


*m^V = 1 WcosmW 


7. 3 A Superposition of Huygens' Sources 

Consider two lines of crossed electric and magnetic current sources 
located on the surface of a dielectric sphere as shown in Fig, 7-3. 

A 

In the x-z plane, there is a filament of +0 directed electric current, 

J (R*) = X(6') 5(R, " - r [6(^)-6(^-,r)]0 , 
e asm©' L -l 

together with a -y directed magnetic dipole strip. 
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FIG. 7-3: A SUPERPOSITION OF HUYGENS' SOURCES ON THE 
SURFACE OF A DIELECTRIC SPHERE. 


J (R 1 ) = -r}„X(0') 
m 2 


6(R* - a) 
asin0' 


+ 6(^' -7r)jy , 


where X(0') is the electric current distribution in the ^ = 0 half-plane. Note 

that at any value of 0', J = -n 0 J > and hence there is a continuous line of 

m 2 e 

approximate Huygens' sources. 

In a similar manner, in the y-z plane the line of Huygens' sources 
consists of an x directed electric dipole strip. 


6(R'-a) 

asin0' 


IT J6(^' - -| ) +6(^' + -| )] x , 


together with a filament of to directed magnetic current 


5 <R') =r? 9 Y(0') 
m 2 


6(R'- a) 
a sin 6* 



Here Y(0') represents the electric current distribution in the <f> = jt/ 2 half- 
plane. 
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Proceeding in a manner similar to that used in Sec. 7.2, it is found 
that the far-zone field radiated by this superposition of sources may be expressed 
by 


E(R) 


ik lR ^ fa . m) , 

7rR e n(n+l) 1 (n + m)! 


- ( A<1>+ W§=T X >+Y* (0 J1 m 

V n Ve /L mn x mn Y_1 omn 
h. ' T r 

f ° (1) \ "1 

+i (b (1) + (6)+Y (0 )] n i , 

yn \/ e r / L. mn x mn yJ emn [ 

where 


X 
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mn y 
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Y' <8 ) 
mn y 


' 0 x P m (cos 0') 

m X(0')d0’ 


sin0’ 


'0 


\0 m+3 m. ... 

Y — r— P (cos 0*) 

(-1) m n . ■— Y(0*)d0* 
sin0' 


^x 9P m (cos 0’) 
n 


90 1 


X(0’)d0’ 


'0 


N0 y ^ 9P m (cos 0') 

(-D — Y(0')d0' 


90’ 


/0 


(7. 17) 


7.4 Numerical Results 

The dipole patterns shown in Figs. 7-4 and 7-5 were computed using 
(7. 11) with the sinusoidal current distribution 


1(0') =sink o a(0 o -0') . 


(7. 18) 
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The upper patterns in these figures show the effect of bending the dipoles into 
arcs with no sphere present. In Fig. 7-4a, the dipole is in the form of a semi- 
circle, i.e., = tt/ 2, and similarly in Fig. 7-5a, the dipole subtends a half- 

angle of 7r/l0 radians. It is noted that these patterns are the same in both the 
forward and back directions and that arcing the dipole tends to reduce the nulls 
in the E-plane patterns occurring at 0 = 7r/2. 

The patterns of Figs. 7 -4b and 7- 5b show the effect of adding a 

dielectric sphere with e = 3. 2. Here the sphere tends to focus the energy into 

r 

the forward direction. 

The current distribution of (7 . 18) was also used in calculating the 
patterns of Fig. 7-6. They illustrate that directional patterns can be produced 
with a small sphere constructed of a material with a relatively high (e^ =7. 5) 
dielectric constant. 

The patterns shown in Fig. 7-7 were calculated using (7. 17) in 2-degree 
intervals with sinusoidal current distributions in both planes. The crossed 
sources each subtended half-angles of 7r/ 6 radians corresponding to a 3-inch 

sphere against the aperture of a C-band circular waveguide. In a comparison of 
Figs. 7-7 and 6-1, it is seen that extending the source had very little effect on 
the resulting patterns. 

Using (7. 17), the patterns of Fig. 7-8 were calculated with a) a point 
Huygens’ source and b) an extended crossed source with sinusoidal currents, 
each subtending half-angles of 7r/ 10 radians. This would correspond to a 5-inch 
sphere operating at 5. 6 GHz. Here it is noted that with the extended source the 
sidelobes are slightly reduced. 



RELATIVE POWER (dB) RELATIVE POWER (dB) 
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FIG. 7-4: RADIATION PATTERNS FOR CURVED HALF- 
WAVE DIPOLE ON SPHERE SURFACE (D = X /tt, 
tan 5 = 0.0). ° 


RELATIVE POWER (dB) RELATIVE POWER (dB) 
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FIG. 7-5: RADIATION PATTERNS FOR CURVED HALF- 
WAVE DIPOLE ON SPHERE SURFACE (D =5\ /? r, 
tan 6 =0.0). ° 



FIG. 7-6: RADIATION PATTERNS FOR CURVED HALF-WAVE DIPOLE 

ON SPHERE SURFACE (D = AX , e = 7.5, tan 6 = 0.0). 

o r 

FIG. 7-7: RADIATION PATTERNS FOR EXTENDED SOURCE (D = 1.27X , 
€^=2.57, tan 6 =0.0, = tt/ 6). ° 






Chapter VTH 

INTERPRETATION BASED ON GEOMETRICAL OPTICS 

8. 1 Introduction 

Physical insight into the way in which the homogeneous dielectric 
sphere focusses the energy from a point source to form the relatively direc- 
tive antenna patterns shown in Chapter 6, may be obtained using the ray-tracing 
techniques of geometrical optics. With this method, it is assumed that the wave- 
length is infinitesimally small, i. e. , ka— ► oo, and the paths of the electromag- 
netic wave propagation are determined using elementary geometry. 

In this chapter, the trajectory of ray paths leaving the sphere are 
determined for various dielectric constants. Next, the amplitude and phase 
across an equivalent plane through the center of the sphere is examined and used 
to explain some of the results of Chapter 6. Finally, we look at the reduction of 
aperture radius caused by the critical angle as a function of source position and 
dielectric constant. 

8.2 Ray Paths 

Figure 8-1 shows the cross section of a homogeneous dielectric sphere 
with an index of refraction n, where n = . Consider a typical ray path 

ABF eminating from a point source located at A; the ray originates at the 
angle 0. and travels undeviated to B. Here the ray is refracted by the air- 
dielectric interface with the transmitted portion leaving the sphere at an angle 
0' measured from the horizontal. 

The plane of incidence is defined as the plane containing the incident 
ray and the normal to the surface. Since for this case the surface normal is the 
radius vector, the plane of the paper is the plane of incidence. It follows, there- 
fore, from the law of refraction (Snell's law) that the refracted ray BF must 
remain in the plane of the paper. To find the angle 0', Snell's law is applied at 
the local point B, 

sin0 = nsin0. (8, 1) 

t X 
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FIG. 8-1: GEOMETRY FOR RAY TRACING CALCULATIONS 
and hence 

0' = 20.- sin 1 (nsin0.) . (8.2) 

Since rays are travelling from a dielectric into a region of lower per- 
mittivity, there is a critical angle 0 given by 

c 

0 = sin *(l/n) . (8.3) 

c 

Rays incident at this angle, i. e., 0. =0 , leave the sphere tangent to the surface. 

1 c 

If 0^ > 0^, the rays experience total internal reflection and will not be considered 
here. 

Ideally one would like the dielectric sphere to transform the radiation 
from the point source into a plane wave. This would correspond to 0' = 0 for all 
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values of 0.. For small 6. (8.2) becomes 
i l 


0' = (2-n)0. 

l 


(8. 4) 


and hence the emergent rays for this case are parallel to the axis when n = 2. 

j 

This value of index is in agreement with Born and Wolf (1970, p. 162), since for ' 
a spherical lens of diameter D the paraxial focus is given by the formula 


D n 
4 n-1 


(8. 5) 


and therefore lies on the surface when n * 2 corresponding to a dielectric 
constant of 4. 

The deviation angle 0' as given by (8. 2) is plotted in Fig. 8-2 as a 
function of incident angle for various dielectric constants. It is observed that 
except when 0. = 0, the deviation angle is strictly positive or strictly negative for 
for e p < 2 or €. > 4, respectively. It appears, therefore, that the sphere will 
have its maximum focussing effect when the dielectric constant is somewhere 
between these two values. Note also from Fig. 8-2 that the critical angle de- 
creases with increasing dielectric constant. Figure 8-3 shows typical ray paths 
for a sphere with a dielectric constant of 3.0. 


8. 3 Amplitude and Phase Across an Equivalent Plane 

Now an expression is derived for the amplitude and phase across an 
equivalent aperture plane, located at the center of the sphere, in a manner 
similar to that outlined by Bekefiand Farnell (1956). 

In this analysis we consider incident rays in the range 


0 < 0 . < 0 , 

— i — c 


e > 2 
r — 


0 < 0. < tt/ 4 , e < 2 

- i _ / r 


since the critical angle occurs at 7r/4 for = 2. Looking back at Fig. 8-1, 
a typical ray travels from A to B within the dielectric, and then leaves the 
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FIG. 8-2: ANGULAR DEVIATION OF RAYS FROM HORIZONTAL 
UPON EXIT FROM SPHERE. 
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FIG. 8-3: RAY PATHS THROUGH DIELECTRIC SPHERE (e = 3.0). 

r 


sphere at the angle 0'. Reflection from the surface is being neglected. Now the 
sphere is removed and the ray BF is traced back as if it originated from point 
H on the plane PP\ 

The various angles and distances in Fig. 8-1 may be found as follows: 

0 = sin ^(nsin0.) 
t l 

0 ' = 20 .- 0 
l t 

AB - Dcos 0. 

i 

OE = CB = 7T cos 20. 

" 1 

EB = OC = ~ sin 20. 

ct X 


BH = CB sec 0* 
CH = BH sin0' 


( 8 . 6 ) 
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Since the ray ABF suffers no phase reversal at B, the phase angle ^ 
of the ray BF transformed back to H is given by 


i> = [ E 


nAB - BH -D(n 


-$»] 


radians 


(8.7) 


where D is the sphere diameter, k^ is the free space propagation constant, 
and the factor kgD(n - j ) is used to normalize the phase to zero at the center 
of the sphere. Making use of (8. 6) we obtain the desired phase expression 


f~ , / cos20.\H 

j> = k Q D n(cos0 ; - 1) + - (l- 1 


cos 0. 


l 


( 8 . 8 ) 


It is apparent from Fig. 8-1 that the radius H of the intersection of the 
ray with the aperture plane is given try 

H = OC - CH 


and hence 


H = — (sin 20 -cos 20. tan 0’) 

2 i i 


(8.9) 


The amplitude of the power density across the equivalent aperture plane 
remains to be examined. For the present we neglect the transmission coeffi- 
cient associated with the surface of the sphere and consider two concentric cones 
of semi-angles 0 and 0+d0 as shown in Fig. 8-4. The power radiated by the 



FIG. 8-4: GEOMETRY FOR APERTURE PLANE AMPLITUDE 
CALCULATION. 
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the source into the space enclosed by these cones appears over the equivalent 
plane in the annular ring of radius H and width dH. 

The incremental areas of the annular rings located at unity radius and 
at the equivalent plane respectively are given by _ . 

dS^ = sin0d0d^ 
dS 2 = HdHdj^ . 


If P(0) is the power radiated per unit solid angle and P(H) is the power received 
per unit angle across the equivalent plane, conservation of energy requires 
(Born and Wolf, 1970, p. 116) 


P(H) 

P(6) 



and thus 


P(H) = P (0) 


sin 0 d0 
H dH ' 


( 8 . 10 ) 


Using (8. 9) together with the relations (8. 6) yields 

2sin0. 2 "I' 1 

P(H) = — — — 2 cos 20. +2 sin 20. tan 0* - cos 0. sec 0 T (2 -ncos 0 sec 0 ) P(0). 

HU L. i . 1 1 l t J 


( 8 . 11 ) 


Although only the primary rays are examined in this analysis, the above 
equation must now be modified to account for the amplitude reduction of these 
rays due to the energy reflected from the sphere surface. 

Since the Huygens' source radiates a linearly polarized field, the two 
principal planes are considered where the electric field is either parallel or 
perpendicular to the plane of incidence corresponding to the E or H planes. 
The power transmission coefficients for these two planes are (Stratton, 1941, 
p. 496) 
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sin 26 i sin 2 e^ 

2 " 2 ( 8 - 12 ) 
sin (0. + 0)cos ( 6 .- 6 ) 
it it 


sin 20. sin 20 

T = __i . (8. 13) 

H sin (0. + 0 ) 
i t 


With this modification, the amplitude distribution A(H) across the line of inter- 
section of the plane of incidence and the equivalent plane will thus be given by 


A(H) = |p(H)t] 1 ^ 2 , (8.14) 

where T is either T or T for the plane of incidence parallel to the E or H 

EH 

planes, respectively. 

Comparing (8. 12) with (8. 13) it is seen that the former contains an 
additional cosine factor in the denominator, and hence the amplitude in the E- 
plane will, in general, be greater than the H-plane. Since the two transmission 
coefficients are identical for 0. = 0^, which would correspond to all rays for 
= 1 or the on-axis ray for > 1, the amplitude distribution will be contin- 
uous at the center of the aperture. 

The results of this section are presented in Figs. 8-5 to 8-7 for a sphere 
of unity diameter. In Fig. 8-5, the phase shift (8. 8) is plotted as a function of 
equivalent plane radius (8. 9) for various dielectric constants. Here a dielec- 
tric constant of about 3.0 tends to yield the minimum phase deviation across 
the aperture. Looking back at Fig. 6-9f it is noted that this same value of 
dielectric constant produced the maximum directivity for the larger diameter 
spheres. That this should be the case is apparent from (8. 8), where we see 
that the phase deviation is proportional to sphere diameter, and hence is a 

critical parameter for large spheres. 

2 

With P(0) = (1+cos 0.) in (8. 11), the amplitude distribution across the 

equivalent plane (8. 14) is plotted in Figs. 8-6 and 8-7 for various dielectric 

constants. We note that for e >2, there is a greater tapering of the amplitude 

r 

across the aperture in the H-plane than the E -plane. In the light of this one expects 
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FIG. 8-5; PHASE ACROSS APERTURE FOR VARIOUS DIELECTRIC 
CONSTANTS. 



Normalized Amplitude 
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FIG. 8-6: E-PLANE APERTURE AMPLITUDE FOR VARIOUS 
DIELECTRIC CONSTANTS. 


Normalized Amplitude 
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FIG. 8-7; H-PLANE APERTURE AMPLITUDE FOR VARIOUS 
DIELECTRIC CONSTANTS. 
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the far-field pattern in the H-plane to have a wider main beam and lower side 
lobes than the E-plane (Silver, 1963). In the patterns presented in Figs. 6-2 
to 6-4, this tends to be the case. 

8. 4 Effective Aperture Radius 

It was noted in Section 8. 2 that rays incident at an angle equal to or 
greater than the critical angle experience total internal reflection. We now wish 
to investigate how the effective aperture radius, which is limited by this critical 
angle, varies as a function of the dielectric constant for the source located on the 
surface or inside the sphere. 

Figure 8-8 shows the cross section of a sphere of radius a. Assuming 



FIG. 8-8: GEOMETRY FOR APERTURE RADIUS 
AT CUTOFF. 


the ray AB is incident on the surface at the critical angle 6 , we seek the 

c 

aperture radius d. As was shown in Section 8.2, the critical angle is given by 

s in 0 = 1/e” = 1/n 

c y r 


(8. 15) 
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where n is the index of refraction. Also from the Law of Sines note 


and hence 


sin0 sin0, 

c 1 


sin0, = — . 

1 nb 


(8. 16) 


Since 


e 0 = e+e , 

2 1 c 


— = sin0,cos0 +sin0 cos 0, 
a 1 c c 1 


(8. 17) 


Using (8. 15) and (8. 16) in the above expression yields 


d 

a 



a I V 

b NTr -1 ; 


(8. 18) 


which is valid for (a/b) < e^. 

Equation (8. 18) is plotted in Fig. 8-9 together with the optimized 
directivity data of Fig. 6-22. Here the effective aperture radius is generally 
about nine-tenths of the sphere radius for the optimized points, and that locating 
the source at the paraxial focus would result in a considerable reduction of the 
of the effective aperture radius. 
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Chapter IX 

CONCLUSIONS AND RECOMMENDATIONS 
9. 1 Conclusions. 

This thesis has been devoted to the study of the radiation from simple 
sources in the presence of a homogeneous dielectric sphere. The dyadic Green’s 
functions of the third and fourth kinds pertaining to this problem have been pre- 
sented and should be applicable to most electromagnetic problems dealing with 
homogeneous spheres. Using the method of scattering superposition, exact 
representations for the electric fields excited by infinitesimal electric or mag- 
netic dipole current sources located inside, outside, or on the surface of a 
dielectric sphere were formulated. 

Because of the simplicity and directional characteristics of the Huygens* 
source, the antenna parameters using this source in the presence of a dielectric 
sphere were studied in detail. Expressions were developed for the radiated 
power, radiation resistance, directivity and antenna efficiency. These results 
may be easily reduced to that for individual electric or magnetic dipoles by let- 
ting the appropriate constants be zero. 

When the calculated radiation patterns for the Huygens* source on the 
sphere surface were compared to measured patterns for plexiglas spheres on 
the end of waveguides , excellent agreement was found. These results confirm 
the findings of Croswell and Chatterjee (1972) that waveguide-excited dielectric 
spheres may be used to radiate patterns which are superior to those from opti- 
mized circular horns of the same cross-sectional area. Extending the point 
Huygens* source into crossed curved lines of sources did not significantly affect 
the radiation patterns. 

Again using the point source, the directivity was plotted as a function of 
sphere diameter for various dielectric constants. It was found that for spheres 
greater than 2. 5 Xq in diameter, a dielectric constant of about 3.0 tended to 
yield the optimum directivity. 
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The dielectric sphere exhibits a resonance phenomenon in which much of 
the energy is radiated into the side lobes with a commensurate decrease in the 
directivity. This effect was most pronounced for the larger, denser spheres. 

It was found that with increasing sphere diameter, the resonance alternates be- 
tween the two principle planes with each new resonance in a particular plane being 
at the next higher mode. The resonance may not be apparent in the radiation pat- 
tern of a lossy sphere, but will be manifest as a reduction in the antenna efficiency. 
This is because the resonance occurs at a high order mode which is quickly 
attenuated by dielectric loss. 

Using the Huygens' source model, it was shown that by appropriately 
adjusting the dielectric constant, radiation patterns with high directivity and low 
side lobes may be produced with the source located inside the sphere. This sug- 
gests that directional patterns may be radiated by a truncated sphere mounted 
on a waveguide, as shown in Fig. 9-1. 

Expressions were developed for the radiation from a curved electric 
dipole in the presence of a dielectric sphere. The bending of the half-wave dipole 
into an arc tended to reduce the sharp nulls in the E -plane patterns. It was 
found that directional patterns may be radiated from a dipole partially wrapped 
around the meridian of a dielectric sphere. This, together with the previous 
results, suggests the use of a monopole mounted on a ground plane covered by a 
dielectric hemisphere as shown in Fig. 9-2. 

Using the methods of geometrical optics, it was shown that a relatively 
small phase deviation occurs across an equivalent aperture plane when the dielec- 
tric constant of the sphere is about 3. This, together with the tapering of the 
field amplitude across the aperture, may be used to explain why the small sphere 
is able to focus the point source into a relatively directive pattern. Since the 
phase deviation is proportional to sphere size, this focussing effect should tend 
to degrade for large spheres. 
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9. 2 Recommendations for Further Research 

There is considerable room for future investigation, both theoretical 
and experimental, of sources in the presence of dielectric spheres. 

According to Luneberg's original theory (1944), an inhomogeneous lens 

2 

with the permittivity profile e (R) = 2 - (R/a) will focus rays from a point 

r 

source on the surface into a collimated beam. It would be interesting to perform 
a detailed numerical analysis of the expressions given by Tai (1956; 1971, Sec. 
46) for the electromagnetic solution of the spherical Luneberg lens. The re- 
sults could then be compared to those presented in this work for the homogene- 
ous sphere. 

Using the method of dyadic Green’s functions, together with scattering 
superposition, exact expressions could be formulated for the radiation from 
stratified spheres. Numerical results could then be easily obtained with a 
computer program similar to that presented in Appendix B. since only the cal- 
culations of the scattering coefficients would need to be modified. Again, the 
data could be compared to that of a homogeneous sphere to see the effect of lay- 
ering on the antenna parameters such as the radiation pattern, resonance, and 
dielectric loss. 

It would be interesting to construct the models shown in Figs. 9-1 and 
9-2 and measure their performance. The data presented in Fig. 6-22 could be 
used to provide the initial design parameters. The antenna of Fig. 9-1 has a 
smaller axial dimension than the full sphere, as well as providing a flat surface 
for mounting the waveguide. 

The radiation patterns of the antenna shown in Fig. 9-2 should be the 
same as that of a dipole within a full sphere because of the imaging properties 
of the conducting ground plane. This structure is similar to that used by Meinke 
(1970), except that here the source is located inside the dielectric hemi- 
sphere. It could be easily mounted on the outside of an aircraft and would not 
i require a radome. 
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As a final remark, antennas are frequently designed without consideration 

of the radome and then modified to reduce the often deleterious effects of this 

< 

protective cap. In view of this work, perhaps the dielectric cover should be con- 
sidered an integral part of the radiating system used to improve the beam- 
forming properties of the antenna. 




FIG. 9-2: PROPOSED ANTENNA FOR EXPERIMENTAL STUDY. 
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Appendix A 

EVALUATION OF [pjjip)]' \(p) - jph Q (p)J' j Q (p) 

Let z^p) and z 2 (p) be two spherical Bessel functions which are 
particular solutions to the differential equation 

p 2 z" + 2pz' +|p 2 - n(n+ l)j z = 0 

Since 

p 2 z" = p(pz)" -2pz' 

(A. 1) may be in the form 

2 


(A. 1) 


p(pz 1 ) ,, +[p - n(n + 1)] z ^ = 0 
p(pz 2 )"+[p 2 -n{n + l)]z 2 =0 


(A. 2) 


(A. 3) 


On multiplying (A. 2) and (A. 3) by z 2 and Zy respectively, and subtracting, 
it is found 


pz 2 (pZ l^ , “ pz l^ pZ 2^" = ° * 


(A. 4) 


We next form the integral 


^ [pz 2 (pz 1 )"-pz 1 (pz 2 )^dp = 0 


(A. 5) 


and integrate by parts yielding 


p2z 2 z 'r p2z i z 2 = c 


(A. 6) 


where c is a constant. 


Letting z = h (p) and z = j (p), where h (p) (p)+iy (p), (A, 6) 
z n in n n n 


becomes 


2 2 

P j n (p)’Mp) -P h n (p) f J n <p) = c . 


(A. 7) 
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T£ *iog the l imit 


as P — >0, 


this result together 


c = -i 


with the relations 

^ W ' “ M'-j 


i n fc) 


' >hnW ' = K fc 3’-\fp) 

in (A. 7) ylelds lhe 


result 




(A. 8) 



Appendix B 

COMPUTER PROGRAM 

Although numerous computer programs were used throughout this 
research, the one of primary significance, together with the necessary sub- 
routines, is listed in this appendix. This program was written in standard 
IBM FORTRAN IV for use with the G level compiler, and was run on an IBM 
System 360/67 computer. 

The MAIN program calculates and prints out the normalized far-zone 
radiation patterns (magnitude, phase, and intensity in dB) in the two principal 
planes for a Huygens' source on the surface or inside a lossy dielectric sphere. 
It also prints the coefficients of the m and n functions, the directivity, and 
provides the pattern data (truncated to -40 dB) in an array suitable for plotting. 
This program may easily be modified to compute data for only the electric or 
magnetic current source, or for the case where the source lies outside the 
dielectric sphere, by making the appropriate changes in the coefficient formulas 
of lines 35 to 39. 

The input parameters to the MAIN program are defined as follows. 


DIAM 

Sphere diameter (X^) 


PER 

Dielectric constant 


TAN 

Dielectric loss tangent (X1000) 


SP 

Distance from source to sphere surface (X ) 

THINC 

Increment in 0 variable (degrees). 



Note: The parameter THINC should divide evenly into 180 degrees. If 

THINC = 1, for example, data is calculated from 6 = 0 to 0 = 180 
in one-degree intervals . 

Subroutine ASLEG is used to generate arrays of associated Legendre 
functions of the first order and their derivatives. The parameters are as fol- 
lows: 

NMAX Maximum value of degree (N) desired 

THINC Increment in 0 variable (degrees). See note above. 
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SLEGEN I 
DLEGEN f 


The names of single precision two-dimensional arrays in 
which the output functions are returned: 


SLEGEN(N, M) 


P^cos M) 
N 

sinM 


DLEGEN(N,M) 


9P* (cos M) 
N 

9M 


where M = (N - 1)THINC . 


Subroutine CBESS generates arrays of spherical Bessel and Hankel 
functions of complex arguments and their derivatives. The parameters are 


as follows. 


BES 

DBES l 
HANK f 
DHANK 


Names of single-precision complex arrays in which the 
output functions are returned: 


BES(N) = j N1 (X) 


DBES(N) = ^ [x iN . [( xi] 

HANK(N) = h^fX) 

N-l 


DHANK(N) = [xhJJ^OO] 


X Double precision complex argument. 

NMAX Maximum value of order (N) desired. 

Double precision arithmetic is used throughout these subroutines, with 
the final output converted to single precision. 



c. 

c 


********** MAIN PROGRAM********* 


OOOl 

000? 

oo r -3 


OOOA 
■000 5 
CO 06 

re 07 

COOP 

00 0 9 

00 1 0 
00 1 1 
001 ? 

0013 

0014 

on c 
0016 
0017 
00 IP 

0019 
00 ?0 
or ?! 
00?? 
Cr?3 
0074 
C ? 5 
00 ? 6 
or? 7 
00 ? 8 
C029 
CC 70 
003 1 


DIMENSION DL EGEN ( 45, 1 8 1 ) ♦ R 1( 18 1 ) , R 2 ( 1 8 1 1 , TH( 1 81 ) » EHM AX ( 2 ) 

OntlRLF PRECISION SLEGEN( 45 » 1 B1 ) 

COMPLEX AMPN( 181) ,BNPN( 1 81 ) , ANDPN ( 1 81 ) t BNDPN ( 181 ) ,EH< 181, 2), 81(45) 

? *DR 1 ( 45 ) , B2( 45) , DB2 ( 45 ) » B3 ( 45 ) , DB 3( 45 ) » HI ( 45 ) , DH1 (45) ,H2(45) , 

30H? ( 45),H3(45),DH3(45) ,AN,BN,CNA, I , AX , AN A , BN A , SER » CER , ON 
COMO (_EX* 1 6 P01 , R02, R03 
1 READ 500,01AM, PER, TAN, SP,THINC 
LK= ? 

X= * OOC 1 *TAN 
PI=3. 14159265358979 
CONV=180./PT 
T = ( 0 . , 1 . ) 

MA=l 80./THINC+1 
0 E R = PE R * ( 1 . ♦ I *X ) 

SER=CSQRT(CER) 

POl = DT AM* PI 
RK 1 = R01 
R0? = SER*ROl 
Rnn=Rn? 

R03r (OIAM-?.*SP)*PI*SER 
C R0I=KIA ,R02=K2A ,R03=K2R 
NM AX - 2 . *R 01+ 6 . 5 
NSUM=NMAX+1 

CALL ASLEG( MSIJM, THIMC , SLEGEN, DLEGEN) 

CALL C8ESS(B1,0BI,H1 ,0H1 , R01 ,NMAX ) 

CALL CBESS( B2,D32,H?,0H?,R02,NMAX) 

CALL CBESS(B3»D93»H3, 0H3 ,R03,NMAX) 

SDAN=0. 

DO 20 j=l,MA 
ANPN ( J ) =0 
RNPN f J ) *0 
ANDPN( J) = 0 
BNOPNf J)*0 
20 CONTINUE 
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O'! %-> 

09 33 
0034 
00 3 5 
0036 

003 7 
0038 
0030 

0040 
00 41 
004? 
00 43 

0^44 

004 5 

,00 46 
0047 
0043 
0049 
09 50 
00 51 
00 5? 
C 9 53 

0054 

0055 
00*6 
00 57 
90 c 3 
00*° 
0060 
00 61 
006? 
00 6? 
on 64 
00 6 5 


P» 1 MT 239 
DO ?? N= 1 ,NMAX 
M=N+ 1 

ANA = -I/<R02*(HUMl*DB2<M)-DHl(M)*B2IM)n 

8 MA=-T/ (R 0 ?*(H 1 ( M)*DB 2 fM)/SER-SER* 0 Hl(M)*B 2 (M)) ) 

CNA=B 3 (M 1 -I*DR 3 (M 1 /R 03 

BN=- RNA*CNA 

AN=ANA*CNA 

C AW AND BN ARE THE CO EFF T F T ENTS OF THE M AND N FUNCTIONS. 
21 CONTINUE 

PRINT 199 , AN, BN 
Y=N 

AX=< ( 2 . 4 Y+ 1 . )/IY 4 (Y + l.m*C 0 .-n**N 
RX = ( 2* Y+l .) /?. 

SDAN=SDAN+RX*(CABSIANM'* 2 +CABStBNl** 2 ) 

C SD*N IS THE DENOMINATOR OF THE DIRECTIVITY EXPRESSION 
DO ?? 1 = 1 ,MA 

ANPNIL) =ANPN I L 1 + AX*SLEGEN ( N, L 1 *AN 
BNPN m=BNPN(L) + AX*SLEGEN(M,L) 4 BN 
ANDPNI L )= AND PM ( L ) +AX*DLEGEN( N, L I* AN 
9 NDPNIL l-BNDPNI L)+AX 40 LEGENIN,L)*BN 
22 CONTINUE 
DO 25 K*.l ,2 
25 FHMA XI K ) = 0 
DO 50 L= 1 ,MA 

EHIL ,1. ) = ANPN ( t ) +BNDPNI L I 
EHCL ,2 ) =- AND PNC L ) -BNPN I L ) 

DO 50 K= 1 , 2 

43 FHMAX(K9=AMAX1(EHMAX(K) , C A8S < F HI L , K I ) ) 

50 CONTINUE 

DTRH= 10 .*ALOG 10 t I CABS I EHI MA, 1 } )** 2 )/SDAN) 

DO 5 5 K = 1 , ? 

DO 55 L=l,MA 

55 EH(L »K 1 = EH( L»K)/FHMAXIK) 

PRINT 302,01 AM, PER, RK 1 ,R 00 ,N,X 
_ PRI NT 305 
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00 66 
,0067 
006 8 
.0069 

0070 

0071 
007 ? 
00 7? 
0074 
00 75 

0076 

0077 

0078 

0079 

0080 
0081 
008? 
00 83 


:C084 

,0085 

i0086 


0087 

0088 
0089 


OQOO 


PR TNT 303, SP 
PRINT 304 ,DI RH,LK 
00 75 L = 1 * MA 
ANGL E=(L-lf*THINC 
TH ( L ) -ANGLE 
AM AS 1=C AB S ( E H( 1. » ! ) ) 

AMAG?=CABS(EH(L,2)) 

ANGl=ATAN2(A!MAG(EH( L , 1 II , RE AL ( EH< L , 1 > ) ) *CONV 
ANG2=ATAN2 ( A IMAG (EH(L,2I ) , REAL ( EH ( L » 2 ) U*CONV 
P1(L)=20.*ALOG10!AMAG1) 

P? ( L )=23.*AL0G13 ( AMAG2) 

PRINT 310 , ANGLE, AMAG1,ANG1 ,R1 (L) , ANGLE, AMAG2, ANG2,R2 ( L) 

IF (Rl(L) .LT.-40.) R 1 ( L ) = - 43 . 

IF ( R2(L) .LT.-40. ) R2(L)=-40. 

75 CONTINUE 

77 PPINT 2 1 3 , EHMAX ( 1 ) , EHMAX ( 2 ) 

ICO CONTINUE 

WRITE! 10, 114) MA, (R1 ( L ) * R2 (I. > »L = 1 , MA) 

C THIS STATEMENT PRINTS ON DEVICE 10 THE NUMBER. OF POINTS TO BE 
C PLOTTED AND THEN THE E- AND H-PLANE INTENSITY ARRAYS. 

GO TO 1 

114 FORM AT (14, (2 4F1 0.4) ) 

30? FORMAT!') ',?1X, 'NORMALIZED FAR ZONE ELECTRIC FIELD FROM SOURCE LOC 
1ATFD ON SURFACE OF DIELECTRIC SPHERE */* 0 ***** SOURCE: HUYGEN 
2 DIAMETER IN WAVELENGTHS: * ,F6. 3, 5X, • DIELECTRIC CONSTANT: » 

3 ,F6. 3/5?X, • AIR KA: ' , E6 . 3 , 11 X , • D I ELECT RT C K A: *,F6.3/ 

4 ?8X, 'NUMBER OF TERMS USED IN SERIES: ',I2,19X 
5, • LOSS TAN: ' ,E9.6) 

303 FORM AT ( 1 7X, * SOURCE DISTANCE FROM SPHERE IN WAVELENGTHS: *,F7.4) 

305 F0RMAT(5?X, » SOURCE I NS TOE SPHERE SURFACE') 

304 FORMAT! //, 6X,'PHI: 0 DEGREES ',46X,'PHI: 90 DEGREES 

5 DIRECTIVITY: *,F6.3,* DB'/'O THETA ELECTRIC FIELD 

7 INTENSITY *,16X,* THETA ELFCTRI C FIELD INTENSITY 

8 »/• (DEGREES) MAGNITUDE ANGLE ! DB ) 

9', 1?X, • (DEGREES) MAGNITUDE ANGLE »,7X,'(DB) 

3 ID F0RMAT(5X,?( F7.2,2X,E10.4, 1X,F7.2,3X,F9. 3,24X) ) 


• ,3X, 12) 
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oc<n 

009? 
00 9? 

0094 

0095 


19R 00PM AT ( 10X,2E12.4,5X, 2012.4, 5X , 2E1 2. 4 , 5X , 2f. 1 2 . 4 1 

200 FORMAT! *0* , 1 4X , * M CQEFFIC 1ENT* ,l' + X t *N COEFF- If. \ ENT • f / 1 5X , ' RE At. « , f>X f 
l * IMAGINARY* , 10X» 'REAL ' t 6X, 'IMAGINARY' ) 

21? FORMAT! 'ONORMAL I ZATIPN CONSTANT: ' »E12.4 t 36X»' NORMAL I Z AT ION CONSTA 
INT: *,E12.4) 

500 FORMAT! 5F10.3» ID 
ENO 



0001 

n 

SUBROUTINE C BESS ( BFS f DBFS *HANK t DH ANK » X ,NMAX ) 


c 

INDEX IS ONE GREATER THAN ORDER OF THE FUNCTION 

•00 ' 2 


COMPLEX BES< IlfDBESn. )*HANK(1)*DHANK(1)» I 

00 0? 


C0MPLEX416 X » JO f J ( 90 1 ♦ P » V ( 7? ) 

0004 


i=n .t l.i 

000? 


MN=NMAX+lfe 

000? 


NQ=NN-l 

0007 


NP=NN-2 

0003 

r 

N0=NMAX+1 

0000 


NR=NMAX+2 

00 10 • 


J(NN1={0. ,0. > 

0011 


j(N0)»u.»i.y 

00 12 


JO=CDSIN( X) /X 

0013 


Vf 11 =-COCOS< Xl/X 

0014 


Y(2)=Y(11 /X-JO 

001? 


DO 10 M=1 » NP 

0016 


N=N3-M 

001.7 

1.0 

J(N1=(?*(N-1 ) +3 ) *J ( N+l l/X-JIN+2) 

00 1 8 


DO 1 2 N=3 ,NR 

00 1 0 

l? 

YfNl =U 7*<N|-1 1-1 )/X)*Y<N-l 5 -YIN-2) 

00 ?0 


p=jo/jm 

00 21 


BF S ( 11=J( 1>*P 

007? 


HANK(i)=?Esm+mm 

00 0 3 


nn 20 L=1 » NO 

■00 ? 4 


K=L+1 

00 7 5 


BES( K1=J( K)*P 

0006 


HANK(K)=BFS(Kl+r*Y(Kl 

0077 


ORES (L)=L*BFS<L)-X*BFSna 

00 7.3 

20 

DHANK(L)=L*HANK(L)-X*HANK( K) 

007 6 


RETURN 

0030 


END 
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II-C * OH KOEEH ‘C5*l-9CA - Ci6I s 13WJO SNllNlHd 1N3 *n» 3AOS ** *fl 


0001 


n oo 2 
00 00 
00 "A 

f)pn«5 
0006 
ooo 7 
00 0 0 
ocoo 
0010 
0011 
00 12 
O' 1 3 
CD 14 

0015 

0016 

0017 

0018 
0010 
00 70 
00 21 
0022 
00 23 

0024 

0025 
0076 
00 27 
0028 
0029 


SUBROUTINE A SLEG ( NMAX, THINC, SLEGEN, DLEGEN) 

C 

DIMENSION DLEGEN 145, 181) 

DOUBLE PRECISION SLEGENI 45 , 1 8 1 ) , X , TH,R AD ,CONV, P ! 

MMAX=1 80. /THING 
MA = M MAX+1 • 

PI=3. 14159265358979 
CONVERT n 80. ‘ 

P,AD= THTNC *CONV 
DO 10 N= 1 ,NMAX 
A=N • 

SL EG ENIN, 1) -A*( A + I. )/2. 

DLEGEN IN, 1)=SLEGEN(N,1) 

OLEGFNIN,MA) =A*( A+l . )/2. *C-1 ) **N 
SLEGENI N,MA) =-OLEGEN(N,MA) 

10 CONTINUE 

DO 20 M=2,MMAX 
L=M-i 

th=l*rad 

X^DCOS(TH) 

SLEGENf 1,M)=1. 

DLEGEN! 1,M)=X 

SLEGENI 2*M)=1 ,5*DSIN( 2.*TH)/DSTN{ TH) 

DLEG EN( 2 » M ) = 3, *DCOS I 2 . *TH ) 

DO 20 N=3,NMAX 

SLESEN ( N , M) ■= U2. *N-1 )*X* SLEGEN (M-l ,M )-N*SLEGEN( N- 2 , M ) ) / ( N- 1 ) 
DLEGEN IN, M)=N*X#SLEGEN(N,M)-(N+1)*SLEGEN(N-1,M) 

20 CONTINUE 
RETURN 
END 
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